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Radiation and Scattering from Printed Antennas on Cylindncally Conformal 

Platforms 

PROJECT SUMMARY 

The goal of this project is to develop suitable methods and software for the analysis of 
antennas on cylindrical coated and uncoated platforms. Specifically, the finite element- 
boundary integral and finite element- ABC methods were employed very successfully for 
this purpose and associated software were developed for the analysis and design of 
wraparound and discrete cavity-backed arrays situated on cylindrical platforms. This is the 
first work which led to the successful implementation of analysis software for such 
antennas, and several significant developments played a major role in this respect. Among 
them, are the efficient implementation of the 3D Green's function for a metallic cylinder, 
the incorporation of the fast Fourier transform in computing the matrix-vector products 
executed in the solver of the finite element-boundary integral system and the development 
of a new absorbing boundary condition for terminating the finite element mesh on 
cylindrical surfaces. Of importance to this project was also the antenna measured data 
provided to us by Mr. Randy Sliva at the Naval Air Warfare Center, China Lake, CA. 

Since there is a dearth of reference data in the literature, the patterns provided by Mr. Sliva 
were essential in validating (to the extent possible) the antenna analysis part of the new 
formulation and associated software. 

As can be realized, the emphasis of our work has so far been on the analysis of uncoated 
cylindrical platforms. Although of interest to this project was also the simulation of coated 
platforms, it was necessary to first complete the analysis of antenna on metallic platforms 
before proceeding with the characterization of patch antennas on the more complex coated 
configurations. Also, it was necessary to generate reference data for such configurations. 
Toward this purpose, a new finite element-boundary integral code is near completion for 
the analysis of flat patches on coated platforms. In this implementation, much effort was 
devoted to developing an efficient algorithm to evaluate the coated platform Green's 
function. The already completed and more versatile finite element-ABC code for 
cylindrically conformal platforms is already capable of simulating patch antennas on coated 
platforms. However, its validation will not be completed until the results of this code (and 
associated formulation) are compared with data based on the reference code and possibly 
measured data to be collected in the coming months. It should be mentioned that the 
implementation of the reference was delayed because of two dimensional reference code 
was first developed as described in the University of Michigan Radiation Laboratory report 
030601 -2-T. 


This report contains three main section and two appendices 
MAIN SECTIONS: 

1. Scattering bv Cavitv-Backed Antennas on a Circular Cylinder 
This section describes the finite element-boundary integral method and its 
implementation for antennas on cylindrical platforms. The Green's function expressions for 
large and small radii of curvature are also presented in this sect ion. H owever, the specific 
details of the resulting system solution in conjunction with the FFT are given in the U-M 
Radiation Laboratory report 030601- 1-T. Of most importance in this section is the 
presentation of scattering patterns for several patch antenna configurations illustrating the 
curvature effects. From these results we conclude that patch arrays on wraparound 

i 



substrates exhibit much higher RCS when compared to the same array of patches 
configured in individual cavities. Also, curves patches are typically of lower Q than flat 
patches of the same size. 

2. Radiation bv Cylindrical Conformal Antennas 

This section addresses the performance of patch radiators on curved platforms. It 
includes comparisons of measured and calculated antenna patterns for patch arrays. These 
results demonstrate the effectiveness of the first ever robust formulation for the analysis of 
non-planar conformal patch antennas. From the results, it is shown that careful placement 
of the radiator patches is required to suppress back lobes in the radiation pattern. An 
important future task is, of course, the characterization of different antenna configurations. 
We can now begin to consider coupling issues and the effect of treatment and curvature on 
the input impedance and antenna RCS. We intend to submit an more detailed report on this 
subject in the very near future. 

3. Finite Element-Boundary Integral Code Manual 

This section serves as a users manual for the subject code. Also a short description 
of the code's operation is included. The effective use of this code should provide the user 
with a good tool for the analysis and design of patch antennas on cylindrical platforms. 


APPENDICES 

Appendix 1 : ... •• • - • 

This appendix gives preliminary results based on the finite element- 
absorbing boundary condition(FE-ABC) formulation code. This code is more versatile 
than the finite element-boundary integral code discussed above. It potentially allows the 
modeling of composite platforms and surface treatments. Also, this formulation, which is 
based on a new class of conformal ABCs, can be generalized to allow modeling of doubly 
conformal platforms. Additional validation and development of this code is necessary 
before it can be transferred to other users. 

A ppendix 2 

This Appendix describes the finite element-boundary integral (FE-BI) formulation 
for antennas recessed in a coated ground plane. The emphasis on a new efficient evaluation 
of the Sommerfeld-type Green's function for the coated grounded plane and its 
incorporation into our existing FE-BI code. The intent is to use this code to generate 
reference data for validating the more versatile FE-ABC code and associated formulation. 
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Circular Cylinder 
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Abstract 

Conformal arrays are popular antennas for aircraft, spacecraft and 
land vehicle platforms due to their inherent low weight and drag prop- 
erties. However, to date there has been a dearth of rigorous analyt- 
ical or numerical solutions to aid the designer. In fact, it has been 
common practice to use limited measurements and planar approxi- 
mations in designing such non-planar antennas. In this paper, we 
extend the finite element-boundary integral method to scattering by 
cavity-backed structures in an infinite, metallic cylinder. In particular, 
we discuss the formulation specifics such as weight functions, dyadic 
Green’s function, implementation details and particular difficulties in- 
herent to cylindrical structures. Special care is taken to ensure that 
the resulting computer program has low memory demand and mini- 
mal computational requirements,. Scattering results are presented and 
validated as much as possible. 
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1 Introduction 


Conformal antenna arrays are attractive for aircraft, spacecraft, and land 
vehicle applications since these antenna systems have low weight, low drag, 
flexibility, and cost advantages over conventional protruding antennas. The 
majority of previous studies pertaining to non-planar conformal antennas 
has been conducted experimentally due to a dearth of rigorous analysis tech- 
niques. Traditional rigorous techniques involve an integral equation and are 
limited in terms of radius of curvature and structural complexity. Some ap- 
proximate methods have been considered but these are restricted in accuracy 
and element shape. 

Recently, the finite element-boundary integral (FEM-BI) method was 
successfully employed for the analysis of large cavity-backed planar arrays 
[1]. The resulting system is sparse due to the local nature of the finite ele- 
ment method, whereas the boundary integral sub-matrix is fully populated. 
However, by resorting to an iterative solver such as the Biconjugate Gradi- 
ent (BiCG) method, the boundary integral sub-system may be cast in circu- 
lant form allowing use of the Fast Fourier Transform (FFT) in performing 
the matrix-vector products. This BiCG-FFT solution scheme ensures 0(N) 
memory demand for the entire FEM-BI system and minimizes the computa- 
tional requirements. 

In this paper, the FEM-BI formulation is extended to scattering by aper- 
ture antennas conformal to a cylindrical metallic surface. In contrast to 
the planar aperture array, the implementation of the cylindrically conformal 
array requires shell-shaped elements rather than bricks, and the required ex- 
ternal Green’s function must satisfy the boundary conditions on the surface 
of the cylinder. In its exact form, this Green’s function is an infinite series 
which imposes unacceptable computational burdens on the method. How- 
ever, for large radius cylinders, a suitable asymptotic Formula is available 
and herein used for an efficient evaluation of the Green’s function. In ad- 
dition, the resulting BI system is again cast in circulant form to ensure an 
C>(N) memory demand and to take advantage of the FFT’s efficiency when 
carrying out the matrix- vector product. 

A primary difficulty in studying cavity-backed antennas mounted on curved 
surfaces is the lack of reference data. In this paper, scattering calculations 
based on the FEM-BI method are compared with data based on different 
techniques. Although such validation is necessarily limited, it provides con- 


2 


fidence in the formulation’s accuracy so that this approach may be used in 
extending the available reference data. 

2 FEM-BI for Circular Cylinders 

In this section, the FEM-BI formulation is developed for cavities recessed 
in an infinite metallic cylinder, having walls which coincide with constant 
p-,4>- or z-surfaces (see figure 1). As usual, the finite element formulation 
permits substantial modeling flexibility, including cavity inhomogeneities, 
lumped loads, super/substrate antenna configurations, or microstrip line and 
so on. 

The FEM-BI approach possesses both low memory and computational 
demand when implemented with a BiCG-FFT solver. Although the system 
of equations associated with the FEM formulation is sparse, the boundary 
integral sub-matrix is fully populated. However, if the aperture mesh is a 
uniform grid, the BiCG-FFT solver may be employed for that portion of 
the system thus retaining (9(N) memory demand for the entire system. In 
addition, the solver require low computational demand since both sparse 
matrix- vector products and discrete convolutions using FFTs require only 
O(NlogN) operations per iteration. 

The FEM-BI formulation begins with the weak form of the vector wave 
equation followed by specification of appropriate vector shape functions and 
dyadic Green’s function. The resulting FEM-BI equations are then used to 
solve for the total electric fields within the cavity and on the aperture (see 
for example Volakis et al. [2]). For the specific configuration at hand, the 
weak form of the wave equation can be written as 

f j V xWi(p,4>,z)-V xWifafrz) 

JVi\ f*r{p,<t>,Z ) 

- kle r (p,<f>,z)Wj(p , <f>,z) ■ Wi(p , 4>,z)^pdpd(f>dz 

+(k 0 a)H a (j)6 a (i) j [ \W i (a,<f>,z)-p(a^,z)x 
J S % Sj 

<J 2 (o, i,z) x • Wj(a > 4>\z‘)\ it'dz'dtdz = ft* + (1) 

In this, Wi are vector basis functions with support over the volume Vi which is 
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Figure 1: Illustration of a typical cavity-backed antenna situated on a metal- 
lic cylinder and the associated coordinate system. 
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associated with the i th degree of freedom, and in a similar fashion, 5, and 
represent aperture surfaces associated with the i th and j th degrees ofJYeedom, 
respectively. The appropriate dyadic Green’s function is denoted by G 2 and it 
has convolutional ( <f> = <f>—<f>, z. — z — z) form when evaluated on the surface 
of the cylinder, p = a. The unprimed coordinates represent the test point 
while the primed ones denote the source point. The free-space propagation 
constant is given by k 0 = where A 0 is the free-space wavelength. The 
cavity is filled with an inhomogeneous material having relative constitutive 
properties e r and p r . The function 8 a (i)6 a (j) is the product of two Kronecker 
delta functions. Hence, it identifies which pairs of unknowns belong to the 
aperture and accordingly contribute to the boundary integral sub-matrix. 

The FEM-BI equation (1) may be rewritten in matrix form as 



[ 0 ] [ 0 ] 
[ 0 ] [ 0 ] 


[E7\ _ {/“'} 
' {0} 


(2) 


where the entries of [.4] are due to the FEM portion of the formulation and 
[Q] is the boundary integral sub-matrix. In (2), E* p and E" 11 denote degrees 
of freedom associated with the aperture and interior fields, respectively. In 
this, ff xt are functions of the external excitation and will be discussed later 
in the paper. 

An important factor in choosing the finite elements for gridding the cav- 
ity is the element’s suitability for satisfying the mathematical requirements 
of the formulation as well as the physical features of the antenna system. 
Traditional node-based finite elements associate the degrees of freedom with 
the nodal fields and have proven unsatisfactory for three-dimensional electro- 
magnetics applications since they do not correctly represent the null space 
of the curl operator and hence spurious modes are generated [3, 4]. In con- 
trast, edge-based elements correctly model the curl operator and therefore 
the electromagnetic fields. In addition, edge-based elements avoid explicit 
specification of the fields at corners where edge conditions may require a sin- 
gularity. Jin and Volakis [6] presented edge-based brick elements which are 
convenient for rectangular-type structures and cavities. For cavities residing 
in a circular cylinder, shell elements are the natural choice. 

Cylindrical shell elements possess both geometrical fidelity and simplicity 
for cylindrical-rectangular cavities. Figure 2 illustrates a typical shell element 
which has eight nodes connected by twelve edges: four edges aligned along 
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Figure 2: Cylindrical shell element. 
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each of the three orthogonal directions of the cylindrical coordinate system. 
Each element is associated with twelve vector shape functions given by 


W u (p,<f>,z) = W 43 (p,<f>,z) = Wp{p.<p,Z\-,<l>i,Zu-) 

w s6 (p, <i>, z) = Wp{p, <t>,z\-,<f> r ,z b , -), w 87 (p, <f>, z) = Wp(p, <p, z-,-,<j>i,z b ,+) 


w u {p, 4>, z) = w^p, <t>, z ; p b , •, z t , +), W 23 (p, <t>, Z ) = H^(p, <f>, z ; p a , •, z t , -) 

W 5S (p,<f>,z) = -), W'srfp, <£, 2:) = W^p, <?!>, z; p a , -, z 6 , +) 


— * 

Wfsfa, = W z (p, <f>, Z] p b , <f> r, 

<£, z) = VE z (p, <£, z; p h , •, -), 


^1 ^rj *1 ) 

<f>, z) = VE z (p, <£, z; p a , <f> h •, + ) 


(3) 


where W)* is associated with the edge which is delimited by local nodes (l,k) 
as shown in figure 2. As seen from (3), three fundamental vector weight 
functions are required for the complete representation of the shell element. 
They are 


Wp(/>, <j>, z ; p, J>, z, s) 
fi^(p,^,z;p,«^,z,s) 
W 2 (p, <f>,z\p, ]>,z,s) 


S£b (<ft - j>)(z - z) . 
ah p 

^\p - p){z - z)4> 
7~{P - p)(<t> - 4>)z 


(4) 


where the element parameters (p 0 , p b , <f>i, <j> r , z b , z t ) are shown in figure 2, t = 
Pb — p a , & = <f> r — <f>i and h = z t — z b . The 4-term which appears in the 
definition of the p-directed weight (4) is essential in satisfying the divergence 
free requirement, i.e. so that V • Wj = 0 1 . Note that as the radius of the 
cylinder becomes large, the curvature of these elements decreases, resulting 
in weight functions which are functionally similar to the bricks presented 

v Wj(p,<j>, z) will only satisfy this requirement within the volume of the element. These 
weighting functions introduce artificial charges on the faces of the element and are not di- 
vergenceless at element interfaces. This is allowable since these elements do not guarantee 
normal field continuity across the element faces. 
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by Jin and Volakis [6]. Having specified the vector basis functions, we may 
proceed to develop the matrix entries for the system (2). 

The FEM-BI system is composed of two parts: a sparse FEM matrix and 
a fully populated BI sub-matrix as shown in (2). The FEM matrix entries 
are represented by 




' /t 1 )*J 

1 $t 


- Ke r iy t 


(2)0 


(5) 


where constant material properties have been assumed within each element. 
The subscripts (i,j) refer to the row and column of the matrix entry and 
correspond to the test and source edges, respectively. The auxiliary functions 



V x W t (p,<f>,z; Zi,Si)pdpd(f>dz 

= / W.ip^'Zipj^j^j^j) ■Wt(p,4>,z;pi,4>i,z i ,s i )pdpd<t>dz {6) 

JV X 


are identically zero unless both test and source edges share at least one el- 
ement in common, resulting in a highly sparse system. Physically, such a 
system is a consequence of the locality property inherent in a partial differ- 
ential equation formulation. In (6), the direction of the source and test edges 
are represented by (s, t) € {p, (j>, s}, respectively. Since the edges of the mesh 
are aligned along three orthogonal directions, only six combinations of (s,t) 
are required for / ^ and only three such combinations for 1^ and all of these 
are evaluated in Appendix A. Since (6) is symmetric with respect to source 
and test edges, the FEM matrix will also be symmetric. 

A lumped impedance post may be included in the formulation by adding a 
term to (1) and equivalently to (5); surface or sub-surface metallization layers 
may also be modeled. Radially oriented lumped loads are approximated 
in the FEM-BI formulation by a filamentary load located at ( <t>L,ZL ) [2]. 
Such posts have length /, cross-sectional area s and impedance Z £,. The 
contribution to [.4] is given by 


Aij = jk 0 Z o -!=- f — — — ^lw i (p,<f>,z)W j (p,<t>,z)pdpd<t>dz{ 7) 
sZl Jv p 

which may be readily evaluated in closed form. In addition, infinitesimally 
thin metallization layers may be represented by simply fixing a priori the 
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weight coefficients to zero for weights associated with edges which are tan- 
gential to the metal. This is a consequence of using a total electric field 
formulation. The symmetry and sparsity of the FEM system [.4] is main- 
tained after the addition of these loads while the BI system [£] remains fully 
populated and symmetric. 

The boundary integral provides an exact boundary condition for mesh 
closure and its construction relies on a cylindrical dyadic Green’s function. 
The entries of the boundary integral sub-matrix are 


Gij = ( k Q a ) 2 / / W t (a,<f>,z\pi,^i,Zi,s,) ■ 

J S, J Sj 

p(d, 4>, z ) x G 2 {a, <}>, z ) x p(a, 4>\ z') 
■W,(a, 4>\ z j pj , <f>j, Zj, Sj)d<f) dzdcj)dz 


( 8 ) 


where the weight functions are given by (4) and evaluated at the surface 
p = a. In (8), the dyadic Green’s function (G2) satisfies both the radiation 
condition and the Neumann boundary condition at p = a. This dyadic 
Green’s function may be expressed exactly [8] 


G**(a,lz) 


G zz (a,<j>,z) = 


G**(a,j>,z) = 

1 r°° 

Wr „£o ■/-» 


_±_ f. /» (t,Y 1 aSM 

( 2 tt ) 2 J - 00 \k 0 ) 7 HT\i) 

_jl v r f n M 


1 

'K m h) 

( nk z N 

\ 2 ^ 2) (7)1 

7 


\k 0 ak PJ 



e Hn*~k,z) dkz 


e iW-krz) dkz 


e j(n*-krz) dkz 


(9) 


where 7 = k p a and k p = — k%. However, for large radius cylinders, (9) 

is computationally prohibitive. In these cases, which are of main concern 
in this paper, it is advantageous to employ an asymptotic expression for G 2 
[9, 10, 11, 12]. These employ a creeping wave series expansion of which only 
the two direct path contributions (see figure 3) are retained. The formula 
due to Pathak and Wang [9] 

G zz {a , 4>, z) ~ -yV--’* 0 * J ( cos 2 9 + q{l - q)( 2 - 3 cos 2 6)) v{(3) 
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G* z (aJ,z) ~ 
G u {a,<j>,z ) ~ 


y^< 7 e ;to * sin 0 cos# | (1 — 3^(1 — q)) r(/?)| 

— ^e - - 7 * 0 *! (sin 2 # + q(l — q)( 2 — 3 sin 2 #)) v(0) 

+q\,ec‘‘e(u(l3) -»(/?))] | 

( 10 ) 


where 0 = ks j 3 and q = ^ has proven quite accurate. In the defi- 

nition of 0, s is the usual geodesic path length (s = yj[a&0 + z 2 ) and # is 
the direction of the trajectory (O = tan~ l ). Depending on which of the 
two direct paths (shown in figure 3) is used, $ = 4> or $ = 2n — 0. The soft 
and hard Fock functions, u(0) and v(0) respectively, are characteristic of 
on-surface creeping wave interactions and have been extensively investigated 
by Logan [13]. Although computation of the Green’s function (10) is now 
tractable, evaluation of (8) must be done so that a discrete convolutional 
system is maintained and the singularity of (8) at s = 0 is properly treated. 

Care must be taken in evaluating (8) so that the overall storage require- 
ment remains 0(N) and the singular integrals of (8) are accurately computed. 
If uniform zoning is used, the resulting sub-matrix ([<?]) is block Toeplitz and 
hence amenable to solution using the BiCG-FFT method. For the non-selfcell 
contributions, mid-point integration may be used while a regularization pro- 
cedure must be employed for the self-cell. Bird [12] noted that (10) recovers 
the metallic screen Green’s function when 0 = 0 within the available ap- 
proximation order. This suggests that (8) may be regularized by adding and 
subtracting from (10) the function 

— - r= v vl e~ jk ° R i 

2 G,(aJ,z)= I + ; R= f-r| (11) 

which is the free-space dyadic Green’s function multiplied by two. The re- 
sulting regularized Green’s function (curvature contribution) is given by 

G zz (a,lz) ~ - J ^qe-> k °*{( c os 2 $ + q(l-q)(2-3cos 2 e))[v(0)-l}} 
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G*-(a,4,,s) 


^~qe ’ ko ’ sin 6 cos (1 — 3^(1— q)) [t'((9) — lj 1 

~ 2* ,C l («n 2 0 + ,(1 - q)(2 - 3 sin 2 #)) [ !,(/?) - 1] 

+, [sec 2 * («(/?) -»(/!))]} 

( 12 ) 


and since it is no longer singular it may be evaluated numerically. The planar 
contribution may be calculated in the manner described previously by Jin 
and Volakis [5]. The FEM-BI matrix has now been fully developed and it 
remains to specify the excitation function for external sources. 


3 Plane Wave Excitation 

Plane wave excitation of the geometry is considered in this section for scatter- 
ing analysis. The use of the exact boundary condition in ( 1 ) allows coupling 
of an exterior excitation field into the cavity. We will describe the form of 
the source functional, f ‘ xi , and discuss its numerical implementation. 

The forcing functional, due to exterior sources (/* r< ) is given by 

fi Xt =jZ 0 k 0 af Wi(a,<f> ,z) ■ p{a,4> ,z) y. H cyl {a,(f>' ,z')d<j> dz (13) 

J Si 

where Wi(p, <j>, z) is the testing weight for the i th row of thermatrix and H^ 1 
represents the magnetic field on the cylinder’s surface in the absence of the 
cavity. A plane wave 


E' 

W 


f,' e -:ko(k'r) 

Y 0 {k' x 


Vo 


p' sin 7 cos 6{ — <f>' cos 7 — z sin 7 sin 0 , 


p jk 0 [p sin 0, cos (<t>— 4> % ) + z cos 0 X ] 


(14) 


is assumed to be incident on the cylinder from the direction where 

7 is the polarization angle and e‘ = 0 'cos 7 + <^‘sin 7 is the electric field 
polarization. The total surface field is given by the sum of the incident 
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and corresponding scattered field from the infinite metallic cylinder [14]. 
Specifically, 


= iHf + iHf (15) 


where 






p jk 0 cos6 t z oo 

- 2 n-j r-j £ 

nk 0 a sin 0, n “ oc 
rc sin 7 cos 0, 


cos 7 


^ ko a sin 0, Hn 2 \k 0 a sin0,) 

j2r o ^^e ifc ° co8 * 1 * V 

nk 0 a 


H^\k 0 a sin 0,) 
+<t>— <t>, ) 


r'(2) 


n=— oo [Hn , (fc o asin0 l )J 


(16) 


is obtained from traditional modal analysis. These expressions may be ap- 
proximated by retaining only a few terms of the series if fc o asin0, is small. 
However, as this parameter becomes large (e.g. for large a and 0, — ► 90°), 
(16) may be replaced with equivalent asymptotic representations similar to 
those considered earlier. Utilizing Watson’s transformation and Fock theory 
[14] in connection with (16), we find that 


Hf 


-Ksiii'/sin 0 , [ J <°>(m® p )]' 

P =1 

2 2 
m l 

j 2 Y 0 COS7- 


gjk 0 cos B{Z ^ g—jk 0 a sin 0 t $p |y(°) ) j * 

# (0) (m$ p ) 


k 0 a sin 0j p=1 

2 

-Y 0 sin 7 cos e,e 3kaCOs8 ‘ z £(-1 )P e -;*o° 


P= 1 


-j -r j 9 (x) (™$ p ) 
k 0 a sin 


( 17 ) 


in which = y - {<f> - <f>i), $ 2 = (<f> ~ <f>i ) - f , m = [ ^ a8 2 mg| J 3 , and 
complex conjugation is denoted by an asterisk. The appropriate far-zone 
Fock functions (<^°\ <7^ and f^) are given by Logan [13]. 
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The asymptotic formulas (17) are quite accurate except in the geometrical 
optics region {<f> & <f> t ). In this case, Goriainov’s [15] expressions 

Hf ~ — K sin a sin cos9,z | e -i*oasm «,*i )j* 

+ e jk oat i n e, cos(*-*,) | G (_ mcos (<£ _ ^.))]- | 


T c Vl 

l 4> 


m- ~ 


j2Y 0 cos o- 


- m ; 2 -- e ifcoco.tf,-x| e -jfc 0 a.ing.» 1 [/(°)( m $ 1 )]‘ 
+e jk oa *in8, [ F (_mcOs{<f>- <£,))]* | 


0 -jk o a sin 8 t $1 


+y o sinacos0,e J * oCOS<, ‘ 2 j 

in<?, c(#-^) J G (_ mcos _ ^)) 


<7 ,0) (™$i) 






(18) 


have been found to be more accurate and can be used instead of (17). The 
Fock functions ( G , G ^ and F) are again defined in Logan [13]. These 
surface field expressions may be used to calculate the entries of the column 
vector efficiently via a numerical evaluation of (13). In particular, 

the modal series (16) is used when k Q a sin 6; < 10 and either (17) or (18) 
for k o asm0i > 10 as appropriate. With the excitation functional and the 
FEM-BI matrix now specified, the BiCG-FFT method [16, 17] may be used 
to determine the unknown electric fields within the cavity. 


4 Scattering 

Once the cavity aperture and volume electric fields have been determined by 
solving (2) for an external excitation, the radar cross section (RCS) may be 
calculated. The far-zone fields may be computed by integrating the aperture 
fields with a suitable Green’s function. In this section we present the relevant 
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formula for calculating the far-zone fields and hence the RCS due to excitation 
by a plane wave (14). 

To determine the far-zone fields, we begin with the integral representation 
for the scattered magnetic field in terms of the aperture fields. We have 

H 3 (r,9,<f>) = jY 0 k 0 a J s G 2 {r,d,<i>-,a,<f>',z') ■ 

[, p(a , <f> , z ) x E(a, <f> , z ') j d4> dz (19) 


with (r,0,<^>) indicating the observation point in spherical coordinates. When 
the observation point is very far from the cylinder, the dyadic Green’s func- 
tion in (19) can be replaced by its far-zone representation 


o-jkoT 


G 2 (r, e , <f>; a, <f>\ z) ~ -j—- [G h 9$ + G 0z 9z + G*+W] (20) 


where the unprimed unit vectors are functions of the observation position 
and the primed ones are functions of the integration point in (19). The 
components of this far-zone Green’s function 


G h ~ 

G Qz 

G H ~ 


j 2k 0 cos 0 jkoCOs9z > ^ n in^+U-4 

(27T ) 2 ( k 0 a sin 0 ) 2 n =-oo Hn 2 \k 0 a sin 6 ) 

J ^ C J kg cos Oz e in(f )) 

(27r) 2 a n=-oo Hn 2 \k Q a sin 0) 

i 9 t 00 1 / 

e jk 0 cos6z y^ e J*(f +(4>-4 > )) 

(2^) 2 a sin 6 H {2 \k 0 a sin 0) 


( 21 ) 


are determined by a mode matching procedure. As one might expect, these 
series converge rather slowly for large k Q a sin 0. They must therefore be recast 
in another form by employing Watson’s transformation and Fock theory as 
was done previously (17). In doing so, we obtain 

G H ~ - COS0 e^ocosez yv_i y e -}k o **m 0 * p L( 0 )( m * ) _ )] " 

47 r l k Q a sin 9 

G 6z ~ - k ° S1U0 e jkoCOB6z g—jk 0 a sin 6$ p [ g (0)( m $ )] * 

4jr p= 1 1 J 

2 2 

G U ~ — e jk 0 cosez' e -jk o asm0<t p VT (22) 

2<z7rsin0 1 1 
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where the Fock functions are the same as those used with (17) due to reci- 
procity. As was the case for the plane wave source, Goriainov’s [15] approx- 
imations are more accurate in the geometrical optics region ( <j> % <f>) and 
similar expressions may be obtained for (22) as was found for (17). The far- 
zone scattered field can be computed numerically by using (19) and either 
the series or asymptotic formula as appropriate. Having done so, the RCS is 
calculated from 




lim 47T r 2 

r— * CO 


|fr(r,M)| 


(23) 


Above we presented a FEM-BI formulation suitable for modeling cavity- 
backed structures embedded in a circular cylinder. Next, we consider a few 
numerical calculations aimed at validating this formulation and in giving 
us an appreciation on how the cylinder’s curvature influences the scattering 
parameters. 


5 Results 

Having solved for the electric fields induced by an incident plane wave, the 
resulting RCS data must be validated with known results. As previously 
mentioned, available measured or computed data is rather scarce and as 
a consequence, we are forced to rely on limiting cases in order to validate 
this work. As the radius of curvature decreases, a cylindrical-rectangular 
cavity will approximate a planar-rectangular cavity. Another limiting case 
involves comparison of an elongated 3-D cavity with a corresponding 2-D 
cavity for normal incidence (0, = 90°). Finally, we may compare our infinite 
cylinder results with a finite Body of Revolution (BOR) model for certain 
polarizations and angles of incidence. We begin with the quasi-planar case. 

The first validation effort for scattering by cavity-backed patch antennas 
relies on the fact that a small patch on a very large radius cylinder is quasi- 
planar and approximates rather well an equal sized planar patch. For our test 
we chose as a reference a planar 1.448 x 1.083 patch residing on a 2.89 x 
2.10" x 0.057" cavity filled with a dielectric having t T = 4. The equivalent 
patch on a 10A cylinder is 6.46° x 1.083” residing on a 12.90° x 2.10" x 0.057” 
cavity. Figure 4 shows the results for the patch on a large radius cylinder 
with corresponding data for the planar cavity-backed patch. Clearly, the two 
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Monostatic RCS (oA ) [dB] 



0.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 80.0 90.0 


Angle (8) [deg] 


Figure 4: Comparison of RCS for a planar patch (1.488" x 1.083") residing on 
a 2.89 x 2.10 x 0.059 cavity filled with e r = 4 dielectric and a corresponding 
quasi-planar patch on a large radius (10A o ) cylinder. 
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RCS patterns are in excellent agreement, and although figure 4 illustrates 
only monostatic scattering in the 4> = 0° plane, additional runs for normally 
incident monostatic scattering and various bistatic situations yield similar 
agreement. 

Comparisons may also be made for elongated cavities and 2-D MoM re- 
sults. Long narrow cavities have very little axial interaction for principal 
plane (9 = 90°) excitation and therefore results based on this formulation 
should compare well with corresponding 2-D data. It is well known that the 
RCS of the 3-D scattering body of length L Ao is related to the corre- 
sponding 2-D scattering of the same cross section via the relation 



Such a comparison is shown in figure 5 for monostatic scattering by a 45° x 
5A x 0.1 A cavity for both principal polarizations. Once again the agreement 
between the two results is excellent, thus providing a partial validation of the 
formulation for highly curved geometries. We remark that similar agreement 
has been observed for bistatic scattering in the 9 — 90° plane. 

The planar approximation eliminates the effects of curvature, which is a 
primary interest in this work, and the 2-D comparisons done above are only 
valid for normal incidence. To consider oblique incidence on a highly curved 
structure, we resort to comparisons with a Body of Revolution (BOR) codej' 
for wraparound cavities. Since the BOR code can only model finite struc- 
tures, we simulate an infinite cylinder by coherently subtracting the far-zone 
fields of the finite structure without a cavity from similar data which includes 
the cavity. Such an procedure mimics common measurement practices and 
was found suitable for near normal incidence and quite acceptable near graz- 
ing incidence in the case of H-polarization (a = 90°). An example calculation 
for the latter case is given in figure 6 where a bistatic scattering pattern is pre- 
sented in the <f> = 0° plane due to a plane wave incident at [9{ = 90°, = 0°). 
Clearly, there is good agreement between the FEM-BI results and data based 
on the BOR formulation. 

The previous comparisons serve to validate the formulation. Having done 
so, it is instructive to examine the effect that curvature has on the scattering 
properties of cavity-backed patch antennas. Consider a 2 cm x 3 cm patch 
residing on a 5.0 cm x 6.0 cm x 0.07874 cm cavity which is filled with 
a dielectric having t r = 2.17. Figures 7 and 8 illustrate the behavior of 
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Bistatic RCS (aA 2 ) [dB] 



Observation Angle (0) [deg] 


Figure 6: Comparison of the RCS computed via the FEM-BI method and a 
BOR code for a 3Ao x 0.1 A 0 air-filled wraparound cavity excited by a normally 
incident H-polarized (a = 90°) plane wave. 
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this geometry as a function of frequency and curvature. Evidently, the 
resonance behavior of this patch is sensitive to curvature for both principal 
polarizations. The frequency response for E-polarization is more sensitive to 
curvature since the radiating surface field component is parallel to the long 
side of the patch and cavity. If the patch and cavity were oriented so that the 
long side is in the (^-direction, the response to H-polarization would exhibit 
greater sensitivity. Such an effect is important to low observable antenna 
designers since they want to operate the antenna in the region of lowest 
RCS. This low return region is a consequence of delicate cancellations due to 
the physical layout of the aperture. Such cancellations are not as complete 
for highly curved structures as they are for planar cavities. 

Conformal antenna designers often use wraparound antennas to achieve 
omnidirectional coverage. Two different configurations are typically used: a 
wraparound cavity where the cavity is filled with a single continuous collar 
of dielectric and discrete cavities symmetrically placed around the circum- 
ference of the cylinder. These two configurations are shown in figure (9). 
Since near resonance, the radiation properties of these two types of antennas 
is identical, any RCS advantage which one might possess could govern the 
appropriate choice of arrays. Figure 10 compares the E-polarized monos- 
tatic scattering at 3 GHz in the 6 = 90° plane for a wraparound cavity and 
four discrete cavities; where the patches and cavities are identical to those 
used in the previous example. Not surprisingly, the wraparound structure 
has a higher return due to coupling within the substrate. However, since in 
this case the scattered field is due to the z component of the surface field 
(^-directed magnetic currents), both cavities yield large scattered fields in 
the four directional lobes. Figure 11 is the corresponding comparison for 
H-polarization. In this case, the scattered field is attributed to the <f> compo- 
nent of the surface fields (z-directed magnetic currents). Therefore, substrate 
modes diffract near the patch resulting in discrete lobes for the discrete array 
while creeping waves shed isotropically for the continuous wraparound cav- 
ity. Low observable designs will favor discrete cavity arrays over wraparound 
cavities since the scattering may be channeled in preferred directions and 
the overall scattering level is consistently lower. A final example is shown in 
figure 12 where we observe that other than the expected higher scattering 
from the wraparound cavity, the scattering behavior of the two arrays is very 
similar. 
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RCS (aA ) [dB] 


Frequency [GHz] 


Figure 7: RCS frequency response for a 2 cm x 3 cm x patch residing in a 
5 cm x 6 cm x 0.07874 cm cavity with c T = 2.17 as a function of curvature 
for E-polarization (a = 0°). 




Frequency [GHz] 


Figure 8: RCS frequency response for a 2 cm x 3 cm x patch residing in a 
5 cm x 6 cm x 0.07874 cm cavity with t T — 2.17 as a function of curvature 
for H-polarization (a = 90°). 






(a) 



Figure 9: Illustration of two types of arrays: (a) wraparound array; (b) 
discrete cavity array 
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Monostatic RCS [dBsm] 



Figure 10: Comparison of E-polarized monostatic RCS at 3 GHz for a four 
patch array placed on a wraparound collar or in four discrete cavities. The 
patches and cavities are identical to the one used in figure 7. The observation 
plane is 8 = 90°. 
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Figure 11: Comparison of H-polarized monostatic RCS at 3 GHz for a four 
patch array placed on a wraparound collar or in four discrete cavities. The 
patches and cavities are identical to the one used in figure 8. The observation 
plane is 6 = 90°. 
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Monostatic RCS [dBsm] 



Angle (0) [deg] 


Figure 12: Comparison of H-polarized monostatic scattering at 3 GHz by a 
four patch array placed on a wraparound collar or in four discrete cavities. 
The patches and cavities are identical to the one used in figure 8. The 
observation plane is <f> = 0°. 
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6 Conclusions 


In this paper, we have presented a Finite Element Method - Boundary Inte- 
gral technique suitable for electromagnetic scattering calculations for cavities 
embedded within a circular, metallic cylinder. This formulation is analogous 
to the FEM-BI approach used by Jin and Volakis [1, 5, 6] and may accord- 
ingly be used for the analysis of scattering by a large array of cavity-backed 
patch antennas. These cavities need not be identical, periodically spaced or 
homogeneously filled and may in fact may possess lumped impedance loads 
or surface metallization layers. The FEM approach employs vector finite el- 
ements which properly represent the electromagnetic fields and possess high 
geometrical fidelity for cylindrical-rectangular cavities. Such elements were 
presented and are analogous to the bricks used for modeling rectangular 
cavities. In addition, we presented an efficient method for evaluating the on- 
surface and far-zone dyadic Green’s functions. The presented formulation is 
amenable to solution using the BiCG-FFT method provided uniform zoning 
is used across the aperture and as a consequence, this implementation has 
low computational and memory demand. We have presented some validation 
of this work with appropriate limiting cases which provides further archival 
reference data. In addition, we showed how this formulation may be used to 
influence conformal antenna designs. 
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Appendices 


A FEM Matrix Entries 

The matrix entries for the FEM portion of the system (5) are given in this 
appendix assuming that the cylindrical shell elements (4) are used in (6). 
These integrals are given by 

iff = -^r ru-iM-itW 

Ct J 4>i 

ifn! = ^^2 (4 {fit ~ fit) + 3 (p* + Pt) {fit ~ fib) + 2^^ {fit ~ fit)') + 

{fib ~ fit) ~ 2t (fi» + Pt) + P>Pt In / ( 2 “ *•)(* ~ Z( ) dz 

= ~^r J\p- h){p- fit)dp 

l( zz - [« Q (p? - pt) - t (ps + Pt) + PzPt In j + 

\ {pI ~ Pi) (<f> ~ 

'ff = 

Iff = ^ [j (pi - Pi) + i (P. + M (pi - Pi) + \p.p, (pi - pi) ] x 

/ (z - z s )(z — z t )dz 

Jib 

Iff = [j (Pi - rf) + 5 (A + Pi) (pI - Pi) + \p.P, (Pi - Pi) ] x 
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( A— 1 ) 


/ {4> - 4>»)(<f> - 4>i)d<t> 

J<t>i 

Each of the above unevaluated integrals is of the form 

ll (< - f.) (f - £.) H = \ (i 2 - V 2 ) ({, + £,) + | ( V 3 - i 3 ) + (U-L) 

(A-2) 

The integrals /,* ’ 2 are used in the assembly of the FEM portion ([*4]) of 
the system. 


B Fock Functions 


The asymptotic form of the dyadic Green’s function with observation both on 
the surface of the cylinder and in the far field involves Fock functions. These 
have been extensively studied and tabulated by Logan [13]. The numerical 
evaluation of these functions are performed either for small arguments or 
large arguments. 

The on-surface Fock functions used in this paper are 


v(0 

«(0 





w ^-itUr 


7T Jcoe-M W2( r ) 

j‘3*/4 _£!_ r ^ll e -Kr dT 
Jooe-J 7 */* wo(r) 


(B-l) 


where ^(r) and its derivative w' 2 {t) denote Airy functions of the Second 
Kind. For small arguments (£ < 0.6), the asymptotic expansion of (B-l) is 
given by 

v (0 ~ + +... 

u(0 ~ 1-0 - -i- + + ... (B-2) 

while a rapidly converging residue series is used for £ > 0.6 



n=l 
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10 

= 2e J * y/ir{2 £ ( r n) _1 e ~^ Tn 

n=l 


(B-3) 


where r n and are zeros of w 2 (r) and w' 2 (t), respectively. Those zeros are 
given in the following table 


Table B-l 

Zeros of the w 2 (r) and u> 2 (t) 
t„ = |r n |e _ and r' n = \ T n \e~ 3 ^ 

n 

T n 

kT 

i 

2.33811 

1.011879 

2 

4.08795 

3.24819 

3 

5.52056 

4.82010 

4 

6.78661 

6.16331 

5 

7.94413 

7.37218 

6 

9.02265 

8.48849 

7 

10.0402 

9.53545 

8 

11.0085 

10.5277 

9 

11.9300 

11.4751 

10 

12.8288 

12.3848 


The far-zone Fock functions are given by 


9 {l) (0 

/ (,) (0 

F"( 0 


j l f e^ T 
y/n Jr w'i(t) 
j 1 r e^ T 

y/H Jr u>i(t) 

g (l) (ty£ 

f(0 (Oe ^ 


dr 

dr 


(B-4) 


where ifi(r) and its derivative u4(t) denote Airy functions of the First Kind 
and the integration contour is given by Logan [13]. These functions, g^{£), 
and /^(£), may be calculated using 
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Z < -1.3 

1.39937 + «) m -1.3 <£<0.5 
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m = l 
#ccr'(m)( 


ml 


r'j a' (m)Ai(m) 


0.5 < £ < 4.0 


= 1.8325e[- ,0 “ 23 - J °' 5094K -^] 


z > 4.0 


9 {1) (0 = -j2.0\Z 2 +j 


2 . 0.25 
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0.25 \ ,d 

e J 3 

i 4 ) 


= £ ^ «r- 


m 


£ < - 2.8 
2.8 < Z < 0.5 


771 — 1 

10 [Ko(m)(] 

0.5<{<4.0 
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(B-5) 


Si ^'(m) 

= -1.8325 (0.8823 - j0.5094 + j£ 2 ) e h a8823 -J°- 5094 )«->V] ^ > 4 g 

(B— 6) 


/ (0) (0 = 


(, 0.25 

( 1 “ + 


0.5\ 

rj 6 3 


e < -i-i 


= 0.77582 + e~ jn/3 ]T «) m - 1.1 < £ < 0.5 


_i m! 


771 = 1 


= e-'"' 3 £ e 


10 [«or(m){] 


= 0.0 Z > 4.0 


0.5 < Z < 4.0 


(B— 7) 


with constant k = e j5,r / 6 and the coefficients for (B-5) and (B-6) given in 
the following table 
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Table B-2 

Constants for (B-5) and (B-6) 

m 

c(m) 


/h(m) 

i 

0.7473831 

1.01879297 

0.5356566 

2 

-0.6862081 

3.2481975 

-0.41901548 

3 

-2.9495325 

4.82009921 

0.38040647 

4 

-3.4827075 

6.16330736 

-0.35790794 

5 

8.9378967 

7.37217726 

0.34230124 

6 

56.1946214 

8.48848673 

-0.33047623 

7 


9.53544905 

0.32102229 

8 


10.52766040 

-0.31318539 

9 


11.47505663 

0.30651729 

10 


12.38478837 

-0.30073083 


The corresponding constants for (B-7) are given as 


Table B-3 

Constants for (B-7) 

m 

c(m) 

a(m) 


i 

1.146730417 

2.33810741 

0.70121082 

2 

0.86284558 

4.08794944 

-0.80311137 

3 

-2.0192636 

5.52055983 

0.86520403 

4 

-9.977776 

6.78670809 

-0.91085074 

5 

-14.59904 

7.94413359 

0.94733571 

6 

49.0751 

9.02265085 

-0.97792281 



10.04017434 

1.00437012 



11.00852430 

-1.02773869 



11.93601556 

1.04872065 



12.82877675 

-1.06779386 
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Abstract 

Conformal antenna arrays are popular antennas for aircraft, space- 
craft and land vehicle platforms due to their inherent low weight, cost 
and drag properties. However, to date there has been a dearth of 
rigorous analytical or numerical solutions to aid the designer. In fact, 
it has been common practice to use limited measurements and planar 
approximations in designing such non-planar antennas. In this paper, 
we extend the finite element-boundary integral method to radiation 
by cavity-backed structures in an infinite, metallic cylinder. The for- 
mulation is used to investigate the effect of cavity size on the radiation 
pattern for typical circumferentially and axially polarized patch anten- 
nas. Curvature effect on the gain, pattern shape and input impedance 
is also studied. Finally, the accuracy of the FE-BI approach for a 
microstrip patch array is demonstrated. 
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1 Introduction 


Modern aircraft and missile deigns seek to utilize conformal antenna arrays 
rather than conventional protruding antennas due to their low weight, low 
drag, low cost and flexibility. Although most useful aircraft surfaces possess 
some curvature, the vast majority of available design information is restricted 
to planar elements. Indeed, the literature is rich with approximate [1], nu- 
merical [2] and experimental [3] design and characterization data for planar 
structures. The most common antenna element is a microstrip patch printed 
on a dielectric coated groundplane. Dielectric coated cylinders have also been 
investigated using approximate [4] and numerical [5] approaches. 

Often, it is desirable to enclose each radiating element within a metal- 
lic cavity to suppress parasitic substrate coupling [6]. Approximate models, 
such as the cavity model, are typically not modified to account for the metal- 
lic sidew'alls of the surrounding cavity since the approximations involved in 
this approach limits its operation to resonant patches. Experience has shown 
that a surrounding cavity does not effect the radiation pattern of a resonant 
patch antenna. However, since integral equation formulations are meant to 
operate at any frequency, these formulations need be modified to account for 
the metallic sidewalls. This is done by partitioning the problem into an inte- 
rior cavity region and an open exterior region and enforcing field continuity 
across the aperture. Such an approach requires a complicated dyadic Green’s 
function for the interior region and as is the case with all integral equation 
formulations, the resulting linear system is associated with a fully populated 
matrix and hence imposes a large 0(N 2 ) memory and computation demand. 
Additionally, most integral equation formulations utilize equivalent surface 
currents and are therefore inappropriate for modeling inhomogeneous sub- 
strates. 

An alternative formulation, utilizing the Finite Element-Boundary Inte- 
gral (FE-BI) method, was proposed by Jin and Volakis [7] which was suitable 
for cavity-backed antennas recessed in a metallic groundplane. As with all 
partial differential equation formulations, this approach is associated with a 
highly sparse system which requires only O(N) storage. Additionally, when 
coupled with a Biconjugate Gradient-Fast Fourier Transform (BiCG-FFT) 
solver, the computational burden is only 0(N log(N)). Since this approach 
is a volume formulation, inhomogeneous substrates may be readily modeled. 
This FE-BI method has been successfully used for scattering and antenna 
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performance analysis in planar platforms. 

Recently, the FE-BI method was extended to cylindrical-rectangular and 
wraparound cavities for scattering calculations [8]. New divergence free, high 
fidelity edge-based elements were presented along with an efficient solution 
strategy which exploited an asymptotic evaluation of the appropriate dyadic 
Green’s function as well as the BiCG-FFT solver. The resulting computer 
code was shown to accurately compute the scattering by planar and highly 
curved elements. This paper investigates the accuracy of this FE-BI method 
formulation for antenna performance analysis. Both radiation pattern and 
input impedance calculations will be compared with known results. The 
effect of curvature on the pattern shape, the resonance behavior and the 
input impedance will be explored. 


2 Formulation 

In this section, the FE-BI formulation appropriate for radiation analysis is 
developed for cavity-backed antennas recessed in an infinite metallic cylinder 
(see figure 1). As usual, the finite element formulation permits substantial 
modeling flexibility, including cavity inhomogeneities, lumped loads and mi- 
crostrip feeding lines. 

The FE-BI formulation begins with the weak form of the vector wave 
equation followed by specification of appropriate vector shape functions and 
dyadic Green’s function. The resulting FE-BI equations are then used to 
solve for the total electric fields within the cavity and on the aperture (see 
for example Volakis et al. [9]). For the specific configuration at hand, the 
weak form of the wave equation can be written as 

r f V x Wj(p, <ft, z) • V x Wj(p, <f>, z) 

Jvi\ 

-kle r (p,<f>,z)Wj(p,<f>,z) • Wi(p, <f>,z)^pdpd<j>dz 

+(k 0 a) 2 8 a (j)6 a (i) f f [w i (a,^z)-p(a,<f>,z)x 
J S i «/ Sj 

&(«,#,*) x ««,#',»') • #,(«>', »)] dtdz'itdz = fr< + fr 1 (i) 

In this, Wi are vector basis functions with support over the volume V{ which is 
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Figure 1: Illustration of a typical cavity-backed antenna situated on a metal- 
lic cylinder and the associated coordinate system. 
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associated with the i th degree of freedom, and in a similar fashion, S, and S } 
represent aperture surfaces associated with the i th and j th degrees of freedom, 
respectively. The appropriate dyadic Green’s function is denoted by Gi and it 
has convolutional (<j> = <f> — <j> , z = z — z ) form when evaluated on the surface 
of the cylinder, p = a. The unprimed coordinates represent the test point 
while the primed ones denote the source point. The free-space propagation 
constant is given by ko = where A 0 is the free-space wavelength. The 
cavity is filled with an inhomogeneous material having relative constitutive 
properties e r and p r . The function <5 a (i)6 a (j) is the product of two Kronecker 
delta functions. Hence, it identifies which pairs of unknowns belong to the 
aperture and accordingly contribute to the boundary integral sub-matrix. 

The FEM-BI equation (1) may be rewritten in matrix form as 



[ 0 ] [ 0 ] 
[ 0 ] [ 0 ] 



{ 0 } 

{frl 


(2) 


where the entries of [.4] are due to the FEM portion of the formulation and 
[£?] is the boundary integral sub-matrix. In (2), Ej P and Ej nt denote degrees 
of freedom associated with the aperture and interior fields, respectively. In 
this, }\ nt are functions of the internal excitation and for this paper a radially 
oriented probe feed is considered. 

The matrix entries, [.4] and [£?], are given in a previous paper [8]. In 
addition, the vector elements, dyadic Green’s function evaluation and far- 
zone field formulae are given and are therefore not repeated here. Hence, we 
need only specify the interior source functional in order to model a radiating 
element. That functional is given for general impressed sources as 




Aft w W.*) 

Mr (p, <t>,z) 


+ jk 0 Z 0 J int (p,<j>,z) 


1 


W{(p, (f>, z)pdpd(f> dz 


(3) 


where J ,nt and M ,nt are the impressed electric or magnetic current densities 
representing the sources. For a radially (p) directed probe feed, the impressed 
monopole current located at ((/>,, z 3 ) is given by 


jtnt 


6((f>- <f> t )(z- z , ) 
pl o 

p 


(4) 
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which results in an excitation function (3) 



if the edge-based elements of [8] are used. 

Having specified the finite element and boundary integral matrices as well 
as the internal excitation for those systems, we use the BiCG method to solve 
for the unknown electric fields throughout the computation domain. The FE 
matrix is highly sparse and hence may be efficiently solved using a sparse 
matrix-vector product. It is also important to note that the matrix-vector 
product associated with the boundary integral can be performed using FFTs. 
Hence, the resulting BiCG-FFT solver is highly efficient without consuming 
excessive memory resources. The electric field may now be used to compute 
antenna parameters such as the gain and the input impedance. 

The radiation pattern is computed by integrating the aperture fields with 
the far-zone dyadic Green’s function given in [8] 

H r (r,6,4>) = jY 0 k 0 a f T? 2 (r, 8 ,<f>;a,<f)\z') • 

J s 

[/5(a, <f> , z ) x E(a, <f > , z )] d<f> dz (6) 

with ( v,0,<f > ) indicating the observation point in spherical coordinates. In the 
far-zone, the electric and magnetic fields are related by 

e; = -z 0 hi 

El = z 0 h; (7) 


which is used to compute the antenna gain 


GdB{0,(f>) = 101og 10 


4 * 


+ I0log,o 


AL 


Z 0 R xn \ 


( 8 ) 


where is the wavelength in centimeters, Ri n is the input resistance which 
is given below and E r is the radiated electric field as r — ► oo. 

In addition to the antenna gain, designers are concerned with the input 
impedance of an antenna for feedline matching purposes. The FEM approach 
allows the calculation of the input impedance of the radiating structure in a 
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rather elegant manner. The input impedance is comprised of two contribu- 
tions [10] 


Z,n = Zp + Zd 


( 9 ) 


where the first term is the probe’s self-impedance ( e.g. the probe’s impedance 
in the absence of the patch) and the second term is the contribution of the 
patch current to the total input impedance. The probe self-impedance ac- 
counts for the finite radius of the probe and hence is omitted when a zero- 
thickness probe is assumed. Ignoring the probe-feed’s self impedance, we 
have [10] 

Z\ n = ~JiJ E(Pifa z ) ' Jl nt {p,<f>, z )pdpd(l)dz ( 10 ) 


where the impressed current is given by (4), V{ refers to the volume elements 
containing the probe-feed, the electric field is the interior field at (p,4>, z) 
and I D is the constant current impressed on the probe. Utilizing (4) and 
assuming the edge- based elements presented in [8], the input impedance (10) 
is evaluated as 


ZL = 


[(*■-*)<-*>] <»> 


which must be summed over the four radial edges of the element which con- 
tains the feed. Having specified the FE-BI system, interior excitation func- 
tional and appropriate antenna parameters such as gain and input impedance, 
the formulation may be used to analyze the radiation characteristics of sev- 
eral interesting configurations. 


3 Results 

The FE-BI formulation presented in [8] may be used for antenna performance 
studies by making the modifications outlined in this paper. This method can 
be used to determine the role of curvature in the radiation pattern and the 
input impedance of a cavity-backed patch antenna. Additionally, designers 
are concerned with the effect of the cavity size might have on antenna perfor- 
mance. In particular, since the cavity is used to suppress parasitic substrate 


43 



modes, varying the cavity size can be used to determine the effect of such 
modes on the radiation pattern. 

Two types of antenna elements are investigated and they are shown in 
figure 2 where each patch is aa° x b in size with a denoting the radius of the 
cylinder. A patch whose radiating side walls are axially oriented is termed 
an axially polarized patch and is fed at <f> 3 = Circumferentially (or az- 
imuthally) polarized patches have radiating walls forming constant z-surfaces 
and are typically fed at z, = |. Observation in the 0 = 90° plane is the E- 
plane for circumferentially polarized patched and the H plane for axially po- 
larized elements. The terminology originates with the cavity model for patch 
antennas. We will now characterize a typical cavity-backed patch antenna. 

Several computed and measured antenna patterns have been published for 
patches printed on a coated cylinder. One such patch, which is 3.5 cm x 3.5 
cm, was used by Sohtell [11] to compare the accuracy of the cavity model [4] to 
a surface current integral equation [5]. The measured data was taken at 2.615 
GHz for a metallic cylinder which was 63.5 cm long and had a radius of a = 
14.95 cm. The cylinder was coated with a 0.3175 cm uniform dielectric having 
relative permittivity of e r = 2.32. Data was taken for —180° < <f> < 180° in 
the 6 = 90° plane corresponding to the E-plane for circumferentially polarized 
elements and the H-plane for axially polarized ones. Figure 3 compares these 
measured patterns with data generated using the FE-BI formulation for an 
identical patch placed within a 360° x 7 cm cavity. This wraparound cavity 
best simulated the measured coated cavity. Note that the H-plane patterns 
are symmetric due to the symmetric placement of the feed, whereas the 
E-plane patterns are not symmetric. The placement of the feed was not 
specified in [11]; however, the agreement for the E-plane pattern shown in 
figure 3 indicates that the position used in the FE-BI model (a<^ 4 = -1 cm) 
is reasonable. The feed was placed at z 3 = -1 cm for the axially polarized 
(H-plane) case. 

In a previous paper [8], discrete cavity arrays were found to have a sig- 
nificantly lower radar cross section (RCS) compared to a wraparound array. 
Thus, the size of the cavity had a significant effect on the scattering prop- 
erties of the array. The two antennas presented by Sohtell [11] were placed 
within cavities which were 7 cm high and approximately 30, 50, 90, 180, 270 
or 360 degrees in angular extent. Figure 4 illustrates that azimuthal cavity 
size has little effect on the radiation pattern for a circumferentially polarized 
element. A similar comparison for the axially polarized patch is shown in 
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(b) 

Figure 2: Illustration of (a) a circumferentially polarized patch element; and 
(b) an axially polarized patch element. The radius of the cylinder is denoted 
by a. 
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Normalized Gain [dB] 



Angle (<J>) [deg] 


Figure 3: Comparison of measured [11] and computed data for a circumferen- 
tially polarized element (E-plane) and an axially polarized element (H-plane). 
The antenna (3.5 cm x 3.5 cm ) was printed on a 14.95 cm cylinder with a 
0.3175 cm coating (e r = 2.32). The probe feed was place at (a <j>,,z s ) = (- 
1.0, 0.0) for the circumferentially polarized patch and at (a^ s , z s ) = (0.0, -1.0) 
for the axially polarized antenna. 
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Normalized Gain [dB] 



Angle (<t>) [deg] 


Figure 4: Effect of cavity size on the E-plane radiation pattern of a circum- 
ferentially polarized patch antenna. 





figure 5. The back lobe of the antenna (near <f> = 180°) is very small for 
cavities less than 180° in extent but increases for larger cavities. For cavities 
which lie on the forward face of the cylinder, the substrate modes diffract off 
the cavity walls; an effect which has little influence on the main lobe of the 
pattern. However, for wraparound cavities and cavities which extend into 
the back side of the cylinder, the substrate modes shed like creeping waves 
giving rise to the back lobe. 

Having established the effect of cavity size on the antenna patterns, it is 
instructive to gauge the effect that curvature has on the resonance behavior 
(or gain) of patch antennas. The two antennas were placed in 14 cm x 14 
cm cavities which were place on cylinders with increasing radius. The fre- 
quency was allowed to vary from 2.4 GHz to 2.7 GHz and the peak radiated 
power was recorded at each frequency. For this paper, the radiated power is 
recognized as the first term of (8). Figure 6 illustrates that the resonance fre- 
quency increases with increasing curvature for a circumferentially polarized 
antenna, but the maximum gain is similar regardless of element curvature. 
Note in the cavity model, the radiating edges for a circumferentially polar- 
ized patch are the azimuthal walls of the cavity (see figure 2) which have 
a constant separation regardless of the cylinder radius. However, the axi- 
ally polarized patch has decreasing resonant gain with increasing curvature 
as shown in figure 7. For this patch, radiation is attributed to the axial 
magnetic walls of the cavity model which have increasing angular separation 
with decreasing curvature. These walls radiate strongly away from the pat- 
tern peak (4> = 0°). Accordingly, the gain of an axially polarized antenna 
decreases with increasing curvature. The radiation pattern of a circumferen- 
tially polarized antenna is largely unaffected by curvature as shown in figure 
8 when excited at a resonant frequency. However, the radiation pattern of 
the axially polarized antenna broadens as the curvature increases and this is 
illustrated in figure 9. Once again, both relationships are readily explained 
by considering the effect that curvature has on the orientation of the cavity 
model radiating walls. 

In addition to the gain and pattern of an antenna, designers require the 
input impedance for matching purposes. For the antenna examined above (in 
a 14 cm x 14 cm cavity), the input impedance was calculated from 2.4 GHz 
to 2.7 GHz for various cylinder radii. Figure 10 illustrates that the input 
impedance of a circumferentially polarized patch antenna is not affected by 
curvature while figure 11 shows that increased curvature reduces the input 
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Figure 5: Effect of cavity size on the H-plane radiation pattern of an axially 
polarized patch antenna. 
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Figure 6: Resonance behavior of a circumferentially polarized patch antenna 
for various cylinder radii. 
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Figure 7: Resonance behavior of an axially polarized patch antenna for var 
ious cylinder radii. 
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Radiated Power [dB] 



Figure 8: Variation of the radiation pattern shape with respect to curvature 
for a circumferentially polarized antenna. 





Angle (<J>) [deg] 


Figure 9: Variation of the radiation pattern shape with respect to curvature 
for an axially polarized antenna. 






Figure 10: Input impedance of a circumferentially polarized patch antenna 
for various cylinder radii. The frequency range was 2.4 GHz to 2.7 GHz and 
the cavity size was 14 cm x 14 cm. 
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Figure 11: Input impedance of an axially polarized patch antenna for various 
cylinder radii. The frequency range was 2.4 GHz to 2.7 GHz and the cavity 
size was 14 cm x 14 cm. 
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impedance of an axially polarized patch. This observation agrees the the 
results reported by Luk et. al. [12]. 

In addition to single patches, the FE-BI formulation may be used to de- 
sign microstrip arrays. Such an approach includes mutual coupling between 
elements which is ignored by the cavity model. Furthermore, the FE-BI for- 
mulation consumes less computational resources than a comparable integral 
equation formulation due to the sparsity of the FE matrix. The H-plane 
pattern of a four element array was measured to gauge the accuracy of the 
FE-BI approach. Each element is 2 cm x 3 cm and placed within a 5 cm x 
6 cm x 0.07S74 cm cavity which is filled with a dielectric having e r = 2.17. 
The cylinder is 91.44 cm long and has a radius of 15.24 cm. The cavities 
are placed symmetrically around the cylinder (e.g. a patch is centered at 
0°, 90°, 180° and 270°). Only the patch centered at 0° was excited while 
the remaining patches were terminated with a 50S7 load. The driving patch 
is axially polarized and the feed is located at z 3 = —0.375 cm. Figure 12 
illustrates the excellent agreement between the FE-BI formulation and the 
measured data. 


4 Conclusions 

In this paper, the FE-BI formulation for cavity-backed antennas was pre- 
sented. The data generated by this code for a common cylindrical-rectangular 
patch antenna compared favorably with measured data. Having validated the 
implementation, the FE-BI method was used to study the radiation prop- 
erties of a circumferentially and axially polarized patch antenna. The az- 
imuthal cavity size was found to have little effect on the circumferentially 
polarized E-plane pattern. However, for the H-plane pattern of an axially 
polarized element, the back lobe is significantly larger for cavities which ex- 
tend from the front side to the back side of the cylinder. A wraparound 
antenna exhibited the largest back lobe implying that this lobe is a result of 
creeping wave shedding. Diffraction due to the creeping wave is suppressed 
for cavities not extending to the back side of the cylinder. The presence of a 
back lobe must be considered when designing low observable, jam-resistant 
antennas or antennas on complex platforms (e.g. an antenna near the rear 
stabilizer). Thus, as was the case for scattering reduction, it is advisable to 
use the smallest cavity possible. 
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Figure 12: H-plane pattern for a four element patch array. Each patch is 
2 cm x 3 cm and are placed symmetrically around the cylinder. Only the 
patch centered at 0° is fed while the other patches are terminated with 50fl 
loads. 
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The effect of curvature on the resonance, radiation pattern shape and in- 
put impedance was studied. Both circumferentially and axially polarized an- 
tennas were considered and was found that the resonant frequency increased 
with increasing curvature for both antennas. However, w r hile the gain of the 
circumferentially polarized patch remained constant, the gain of the axially 
polarized patch decreased with increasing curvature. Such an effect is read- 
ily explained by considering which walls of the cavity model radiate for each 
polarization (see figure 2). The radiation pattern for axially polarized anten- 
nas broadens with increasing curvature while the corresponding patterns for 
circumferentially polarized antennas is unaffected by curvature. The input 
impedance of the circumferentially polarized antenna was not affected by cur- 
vature while the input impedance of the axially polarized antenna decreased 
with increasing curvature. We therefore conclude that axially polarized an- 
tennas are strongly affected by curvature while circumferentially polarized 
antennas are not affected by curvature. 

Since the cavity model does not include mutual coupling and the usual 
integral equation formulations are associated with high storage and com- 
putational demand, the FE-BI formulation is attractive for array analysis. 
The H-plane pattern of a discrete four element wraparound array compared 
favorably with measured data. 


58 



References 


[1] Y.T. Lo, D. Solomon and W.F. Richards, “Theory and experiments on 
microstrip antennas,” IEEE Trans. Antennas and Propagat ., Vol. 27, 
PP . 137-145, 1979. 

[2] D.M. Pozar and S.M. Voda, “A rigorous analysis of a microstripline fed 
patch antenna,” IEEE Antennas Propagat., Vol. 35, pp. 1343-1350, Dec. 
1987. 

[3] J.R. James and G.J. Wilson, “Microstrip antennas and arrays, Pt. 1 - 
Fundamental action and limitations,” IEE J. Microwaves, Optics, and 
Acoustics , Vol. 1, pp. 165-174, 1977. 

[4] J.S. Dahele, R.J. Mitchell, K.M. Luk and K.F. Lee, “Effect of curvature 
on characteristics of rectangular patch antenna,” Electronics Letters , Vol 
23, pp. 748-749, 2 July 1987. 

[5] J. Ashkenazy, S. Shtrikman and D. Treves, “Electric surface current 
model for the analysis of microstrip antennas on cylindrical bodies,” 
IEEE Antennas Propagat., Vol. 33, pp. 295-300, Mar. 1985. 

[6] J. Aberle, “On the use of metallized cavities backing microstrip an- 
tennas”, 1991 IEEE Antennas and Propagat. Soc. Int. Syrup., Vol. 1, 
pp. 60-63, June 1991. 

[7] J-M Jin and J.L. Volakis, “A hybrid finite element method for scattering 
and radiation by microstrip patch antennas and arrays residing in a 
cavity,” IEEE Trans. Antennas and Propagat., Vol. 39, No. 11, pp. 1598- 
1604, Nov. 1991. 

[8] L.C. Kempel and J.L. Volakis, “Scattering by cavity-backed antennas on 
a circular cylinder,” submitted to IEEE Trans. Antennas and Propagat.. 

[9] J.L. Volakis, A. Chatterjee and J. Gong, “A class of hybrid finite element 
methods for electromagnetics: A review,” to appear in J. Electromag- 
netic Waves Appl., 1994. 

[10] T.M. Hashaby, S.M. Ali and J.A. Kong, “Input impedance and radiation 
pattern of cylindrical-rectangular and wraparound microstrip antennas,” 


59 


IEEE Trans. Antennas and Propagat., Vol. 38, No. 5, pp. 722-731, May 
1990. 

[11] E.V. Sohtell, “Microstrip antennas on a cylindrical surface,” in Hand- 
book of microstrip antennas, Ed. J.R. James and P.S. Hall, Peregrinus: 
London, pp. 1227-1255, 1989. 

[12] K-M Luk, K-F Lee and J.S. Dahele, “Analysis of the cylindrical- 
rectangular patch antenna,” IEEE Trans. Antennas and Propagat., Vol. 
37, No. 2, pp. 143-147, Feb. 1989. 


60 



N94-2 

User Manual for EXCALIBUR: A FE-BI 
Numerical Laboratory for Cavity-Backed 
Antennas in a Circular Cylinder 
Version 1.2 

Leo C. Kemp el 
Radiation Laboratory 
University of Michigan 
1301 Beal Ave. 

Ann Arbor, MI 48109-2122 

December 16, 1993 




61 


1 Introduction 

The Finite Element-Boundary Integral (FE-BI) technique has been used to 
analyze the scattering and radiation properties of cavity-backed patch anten- 
nas recessed in a metallic groundplane. A program, CAVITY3D, was written 
and found to yield accurate results for large arrays without the usual high 
memory and computational demand associated with competing formulations. 
Recently, the FE-BI approach was extended to cavity-backed antennas re- 
cessed in an infinite, metallic circular cylinder. EXCALIBUR is a computer 
program written in the Radiation Laboratory of the University of Michigan 
which implements this formulation; 

This user manual will give a brief introduction to EXCALIBUR and some 
hints as to its proper use. As with all computational electromagnetics pro- 
grams (especially finite element programs), skilled use and best performance 
is only obtained through experience. However, we will comment on several 
important aspects of the program such as portability, geometry generation, 
interpretation of results and custom modification. 


2 Formulation 

EXCALIBUR implements the FE-BI formulation for cavity-backed antennas 
recessed in metallic circular cylinders. The formulation imposes some re- 
strictions on the geometries which will be discussed. Principally, use of the 
BiConjugate Gradient-Fast Fourier Transform (BiCG-FFT) iterative solver 
requires uniform zoning on the aperture of the cavity. Thus, each surface 
patch has cylindrical-rectangular shape. The basis functions, dyadic Green’s 
function and associated field formulas are given in a paper concerning scat- 
tering by these structures [1]. Modifications of this formulation for antenna 
analysis was given in another paper [2], 

The uniform zoning requirement causes some difficulty in modeling; how- 
ever, with some practice, these difficulties may be overcome. For example, 
the specification of the patch and cavity size must both be expressed by an 
integer number of edges (hence nodes). Thus, if the cavity is twice the size 
of the patch, one has no problem specifying the patch and the cavity with 
the same uniform grid. However, if the ratio of the patch and cavity sizes are 
not integers, discretization may not be possible. This is often the case with 
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a continuous wraparound cavity. Such a cavity is shown in figure 1 along 
with an example of a discrete wraparound array. If the cavity size and patch 
size are not convenient, you must either change the cavity size as possible or 
change the radius of the cylinder. If the radius is changed slightly, it will not 
effect the electromagnetic properties of the structure, but it may allow uni- 
form discretization. Although the restrictions imposed by the uniform zoning 
requirement seems rather stringent, with practice, an antenna designer will 
find that EXCALIBUR is quite flexible. 

3 Compilation 

The first task in utilizing EXCALIBUR is to compile and link the various 
files of the program. The following files are required to run EXCALIBUR and 
are given on the distribution diskette: excalibur.f dyadic./, matrixGenerate.f 
preProcessor ./, rcs.f \ fft.j , gauss. inc, fft.inc and excalibur.inc. These files 
perform the following functions: 

• excalibur.f. Main program, BiCG solver, matrix building subroutines, 
FE-BI subroutine, impedance insert and various auxiliary subroutines. 

• dyadic.f. Compute dyadic Green’s function terms for admittance ma- 
trix, on-surface and far-field Fock functions and gamma function. 

• matrixGenerate.f. Boundary integral and FE matrix terms. 

• preProcessor. f. Geometry /mesh generator. 

• rcs.f. Radar Cross Section, far-zone dyadic Green’s function and plane 
wave excitation functions. 

• fft.f. Forward and inverse Fast Fourier Transform subroutines. 

• gauss. inc: Numerical integration parameters for gaussian quadrature. 

• fft.inc: Include file for 2-D FFT subroutines. 

• excalibur.inc: Main memory allocation file also contains variable dic- 
tionary. 

The distribution disk also contains the following Fast Fourier Transform op- 
tions: 

• fftCooley.f. Classic Cooley-Tukey Radix-2 algorithm. 

• fftSplit Radix. f. Split- Radix algorithm by Sorenson [3]. 
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(a) 



(b) 

Figure 1: Illustration of two types of arrays: (a) wraparound array 
discrete cavity array 
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• fftCRAY.f: Calls the CRAY library vectorized, multitasked 2-D FFT 
routine. 

• fftCONVEX.f: Calls the veclib vectorized 2-D FFT routine. 

Prior to compilation, the user should copy one of these files into fft.f. The 
optimized routines for the CRAY and CONVEX architectures should be used 
when possible while in general the Split-Radix algorithm should give better 
performance than the Cooley-Tukey version. To enhance their efficiency, the 
Split-Radix and Cooley-Tukey files utilize a decimation-in-frequency forward 
transform and a decimation-in-time inverse transform in order to avoid the 
need to perform bit reversal. 

Another file which is included on the distribution diskette is convert- 
ToASCII.f. To save disk space, the geometry information is stored in a 
binary format by preprocessor. The program convertTo ASCII is included 
to produce a human-readable file (ASCII), The nodes, elements, edges, un- 
knowns and other useful information is provided in a easy to read (although 
disk space consuming) format. All the programs are compiled and linked by 
invoking the UNIX make utility. A Makefile has been provided on the distri- 
bution disk. To date, EXCALIBUR has been successfully compiled, linked 
and run on the following architectures/operating systems: SUN, DEC UNIX, 
HP 9000/7xx, IBM RS/6000, Silicon Graphics IRIS, CRAY and CONVEX. 
Three variable must be set within the Makefile : 

• FF: The Fortran compiler name for the architecture. 

• FOPT: The Fortran compiler options i.e. optimization, precision, etc. 

• LOPT : The name of any libraries required for linking. 

The user should uncomment these variables for the target architectures in 
Makefile. EXCALIBUR is constructed by simply typing make at the com- 
mand line, while the binary-to-ASCII conversion program is constructed by 
entering make convert. EXCALIBUR is invoked by typing excalibur at the 
command line while the conversion program is run with the command con- 
vertTo ASCII. Finally, the directory may be cleaned up of all object and 
executable file by typing make clean. 
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4 Geometry Information 

The binary geometry file created by preprocessor contains all the information 
concerning the physical structure under study except for the placement of any 
probe-feeds or lumped impedance posts. Therefore, it is important that the 
user be aware of the geometry entered into the FE BI portion of the code. As 
previously mentions, the convert To A SCII program creates a human-readable 
file from the machine-readable geometry file. 

The first information provided in the resulting ASCII file is the header, 
which contains the number of nodes, number of edges, number of unknowns, 
etc. and an example of the header is shown in figure 2. The next field contains 
the node information. The information given is as follows (see figure 3): 


Each node is associated with a physical location (p,^,z) and a grid location 
(layer, row, column). 

Grid points must be used in the discretization of a geometry since the 
BiCG-FFT solver requires that each node pair lie an integer number of units 
apart. Thus, the distance between two nodes (primed and unprimed) on the 
surface of the cylinder is given by 

R(n,m;n ,m ) = \J{n — n')aA^ + (m — m')Az (1) 

EXCALIBUR distinguishes between grid points and nodes. A grid point can 
be thought of as the intersection of two lines of a piece of graph paper which 
is placed on the surface of the cylinder. A node is a grid point which lies 
within a cavity. The row and column number associated with a node actually 
is the row and column number of the grid point which formed the node. The 
first grid point which corresponds to the lower-left corner of the grid has row 


Column 1: 
Column 2: 
Column 3: 
Column 3: 
Column 4: 
Column 5: 
Column 6: 


Node number 

Radial ( p ) coordinate in centimeters. 

Angular (<f>) coordinate in degrees. 

Axial (z) coordinate in centimeters. 

Layer number from top of the cylinder (aperture). 
Row number from lowest axial coordinate. 


Column number from smallest azimuthal coordinate. 
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Binary filename : test . exc 


ASCII filename : test . ascii 


NODE statistics: 

Total number of nodes: 50 

Number of nodes on the surface: 25 

Number of nodes along the metallic walls: 41 

Number of nodes on surface metallic patches: 0 

Number of nodes which are resistive: 0 


EDGE statistics: 

Total number of edges: 105 

Interior edges: 9 

Aperture edges: 24 

a) substrate edges: 24 

b) resistive edges: 0 

Metal edges (NOT unknowns) : 72 

II II I I I I I UNKNOWNS I I II I I I I I— > 33 


ELEMENT statistics: 

Total number of elements: 16 

Surface Elements: 16 


Figure 2: Geometry header. 
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Node 

rho (cm) 

phi (deg) 

z (cm) 

layer 

row 

column 

1 

1.00000 

-5.00000 

-0.50000 

0 

0 

0 

2 

1.00000 

-2.50000 

-0.50000 

0 

0 

1 

3 

1.00000 

0.00000 

-0.50000 

0 

0 

2 

4 

1.00000 

2.50000 

-0.50000 

0 

0 

3 

5 

1.00000 

5.00000 

-0.50000 

0 

0 

4 

6 

1.00000 

-5.00000 

-0.25000 

0 

1 

0 

7 

1.00000 

-2.50000 

-0.25000 

0 

1 

1 

8 

1.00000 

0.00000 

-0.25000 

0 

1 

2 

9 

1.00000 

2.50000 

-0.25000 

0 

1 

3 

10 

1.00000 

5.00000 

-0.25000 

0 

1 

4 

11 

1.00000 

-5.00000 

0.00000 

0 

2 

0 

12 

1.00000 

-2.50000 

0.00000 

0 

2 

1 

13 

1 . 00000 

0.00000 

0.00000 

0 

2 

2 

14 

1.00000 

2.50000 

0.00000 

0 

2 

3 

15 

1,00000 

5.00000 

0.00000 

0 

2 

4 

16 

1.00000 

-5.00000 

0.25000 

0 

3 

0 

17 

1.00000 

- 2.50000 

0.25000 

0 

3 

1 

18 

1.00000 

0.00000 

0.25000 

0 

3 

2 

19 

1.00000 

2.50000 

0.25000 

0 

3 

3 

20 

1.00000 

5.00000 

0.25000 

0 

3 

4 

21 

1.00000 

-5.00000 

0.50000 

0 

4 

0 

22 

1.00000 

-2.50000 

0.50000 

0 

4 

1 

23 

1.00000 

0.00000 

0.50000 

0 

4 

2 

24 

1.00000 

2.50000 

0.50000 

0 

4 

3 

25 

1.00000 

5.00000 

0.50000 

0 

4 

4 

26 

0.90000 

-5.00000 

-0.50000 

1 

0 

0 

27 

0.90000 

-2.50000 

-0.50000 

1 

0 

1 

28 

0.90000 

0.00000 

-0.50000 

1 

0 

2 

29 

0.90000 

2.50000 

-0.50000 

1 

0 

3 

30 

0.90000 

5.00000 

-0.50000 

1 

0 

4 


<< Remainder of nodes truncated >> 


Figure 3: Node information. 
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= 0 and column = 0. For a wraparound grid, the first grid point is physically 
located at <j> = —180° and once again has row = 0 and column = 0. 

The next set of information provided is the edges which form the cavities. 
The information given is as follows (see figure 4) 






Column 1: 
Column 2: 
Column 3: 
Column 3: 
Column 4: 
Column 5: 


Edge number. 

Left (lower) node forming the edge. 

Right (upper) node forming edge. 

Unknown number (zero indicates a fixed edge (e.g. metal)). 
Orientation (p-, <f>- or z-directed). 

Type of edge (metal, substrate, resistive or interior). 


Each edge is associated with two nodes and hence has an orientation in the 
cylindrical coordinate system. If an edge is metal, since EXCALIBUR uses a 
total field formulation, that edge’s weight is fixed at zero. All other edges are 
unknowns which must be solved using the BiCG-FFT solver. A substrate or 
resistive edge is associated with the boundary integral while interior edges 
contribute only to the FE portion of the system. Currently, resistive cards 
are not implemented in EXCALIBUR. 

The next set of information related the unknowns on the aperture of the 
cavities to their edge number. It also includes the row and column number of 
that edge in the discretization grid. Although this information is useful for 
understanding the mechanics of the BiCG-FFT solver, it is of little interest 
to the general user. The given information is (see figure 5) 

• Column 1: Unknown number. 

• Column 2: Associated edge number. 

• Column 3: Row of this edge in the discretization. 

• Column 3: Column of this edge in the discretization. 


The edges which form each element of the mesh are given next. Each 
cylindrical shell element consists of eight nodes which form twelve edges. 
This information is useful in visualizing the mesh and could be hooked into 
a graphics package to generate a 3-D picture of the mesh. The prototype 
element is shown in figure 6 which displays the node numbering scheme. The 
information given by convertToASCII is (see figure 7) 

• Row 1, Column 1: Element number. 
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Edge 



Node 1 

Node 2 

Unknown 

Orientation 

Type 

1 

26 

1 

0 

rho-directed 

metal 

2 

27 

2 

0 

rho-directed 

metal 

3 

31 

6 

0 

rho-directed 

metal 

4 

32 

7 

25 

rho-directed 

interior 

5 

26 

27 

0 

phi-directed 

metal 

6 

1 

2 

0 

phi-directed 

metal 

7 

31 

32 

0 

phi-directed 

metal 

8 

6 

7 

1 

phi-directed 

substrate 

9 

26 

31 

0 

z-directed 

metal 

10 

1 

6 

0 

z-directed 

metal 

11 

27 

32 

0 

z-directed 

metal 

12 

2 

7 

13 

z-directed 

substrate 

13 

28 

3 

0 

rho-directed 

metal 

14 

33 

8 

26 

rho-directed 

interior 

15 

27 

28 

0 

phi-directed 

metal 

16 

2 

3 

0 

phi-directed 

metal 

17 

32 

33 

0 

phi-directed 

metal 

18 

7 

8 

2 

phi-directed 

substrate 

19 

28 

33 

0 

z-directed 

metal 

20 

3 

8 

14 

z-directed 

substrate 

21 

29 

4 

0 

rho-di rected 

metal 

22 

34 

9 

27 

rho-directed 

interior 

23 

28 

29 

0 

phi-directed 

metal 

24 

3 

4 

0 

phi-directed 

metal 

25 

33 

34 

0 

phi-directed 

metal 

26 

8 

9 

3 

phi-directed 

substrate 

27 

29 

34 

0 

z-directed 

metal 

28 

4 

9 

15 

z-directed 

substrate 

29 

30 

5 

0 

rho-directed 

metal 

30 

35 

10 

0 

rho-directed 

metal 


<< Remaining edges truncated >> 


Figure 4: Edge information. 
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iknown 

Edge 

Row 

Column 

1 

8 

2 

1 

2 

18 

2 

3 

3 

26 

2 

5 

4 

34 

2 

7 

5 

40 

4 

1 

6 

47 

4 

3 

7 

52 

4 

5 

8 

57 

4 

7 

9 

63 

6 

1 

10 

70 

6 

3 

11 

75 

6 

5 

12 

80 

6 

7 

13 

12 

1 

2 

14 

20 

1 

4 

15 

28 

1 

6 

16 

44 

3 

2 

17 

49 

3 

4 

18 

54 

3 

6 

19 

67 

5 

2 

20 

72 

5 

4 

21 

77 

5 

6 

22 

90 

7 

2 

23 

95 

7 

4 

24 

100 

7 

6 


Figure 5: Relationship between unknown number and edge number on aper- 
ture. 
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X 


Figure 6: Cylindrical shell element. 
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• Row 1, Column 2-5: /^-directed edges. 

• Row 2, Column 2-5: ^-directed edges. 

• Row 3, Column 2-5: z-directed edges. 

The final set of information provided is the element parameters as shown 
in figure 6 which includes (see figure 8) 

• Row 1: Element number. 

• Row 2: p a , pb and t — Pb — Pa- 

• Row 3: (f>i , (f> T and a = <f> r — 4>i- 

• Row 3: Zb, z t and h = z t — Zb. 


5 Geometry Generation 

Having reviewed the geometry information provided by preprocessor through 
convertTo ASCII, we are prepared to generate some example geometries. 
Specifically, we shall look at radiation and scattering by a 2 cm x 3 cm 
patch antenna which is placed in a 5 cm x 6 cm cavity, a four element discrete 
array of such cavities and the same radiating array placed in a continuous 
wraparound cavity. These three variations of the same geometry exhibit the 
main classes of structures encountered in practice. Users will find it very 
helpful to check the entered geometry via convertTo ASCII prior to running 
the solver part of EXCALIBUR. In particular, it is useful to discretize the 
cavity without patches present and retain the node information since it will 
be necessary to specify the row and column of the lower-left corner of each 
patch as well as the number of edges along each side. 

EXCALIBUR has a preprocessor module which generates the required 
mesh. It first generates the surface nodes which are then used to create 
the 3-D mesh. An in-house modeling package, such as SDRC IDEAS, may 
be used to create such a surface grid so long as the nodes are an integer 
number of units apart. It would be quite easy to interface such a package 
with EXCALIBUR by replacing the subroutine simpleMesh with a universal 
file reader. However, we have found that the custom mesh routine provided 
with EXCALIBUR ( simpleMesh ) is sufficient for most modeling tasks while 
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Element Edges 


1 

1 

2 

3 

4 


5 

6 

7 

8 


9 

10 

11 

12 

2 

2 

13 

4 

14 


15 

16 

17 

18 


11 

12 

19 

20 

3 

13 

21 

14 

22 


23 

24 

25 

26 


19 

20 

27 

28 

4 

21 

29 

22 

30 


31 

32 

33 

34 


27 

28 

35 

36 

5 

3 

4 

37 

38 


7 

8 

39 

40 


41 

42 

43 

44 


« Remaining 

elements 

truncated » 



Figure 7: Edges associated with each element. 
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Element size parameters: 


Element 

number : 


1 





rhoA 

= 

0.90000 

rhoB 

= 

1 . 00000 

t = 


0.10000 

phiL 

= 

-5.00000 

phiR 

= 

-2.50000 

alpha 

= 

2.50000 

zB 

= 

-0.50000 

zT 

= 

-0.25000 

h 

= 

0.25000 

Element 

number : 


2 





rhoA 

= 

0.90000 

rhoB 

= 

1.00000 

t = 


0.10000 

phiL 

= 

-2.50000 

phiR 

= 

0.00000 

alpha 

= 

2.50000 

zB 


-0.50000 

zT 

= 

-0.25000 

h 

= 

0.25000 

Element 

number : 


3 





rhoA 

= 

0.90000 

rhoB 

= 

1.00000 

t = 


0.10000 

phiL 

= 

0.00000 

phiR 

= 

2.50000 

alpha 

= 

2.50000 

zB 

= 

-0.50000 

zT 


-0.25000 

h 

— 

0.25000 

Element 

number : 


4 





rhoA 

= 

0.90000 

rhoB 

— 

1.00000 

t = 


0.10000 

phiL 

= 

2.50000 

phiR 

= 

5.00000 

alpha 

— 

2.50000 

zB 

= 

-0.50000 

zT 

= 

-0.25000 

h 

— 

0.25000 

Element 

number : 


5 





rhoA 

= 

0.90000 

rhoB 

= 

1.00000 

t = 


0.10000 

phiL 

= 

-5.00000 

phiR 

= 

-2.50000 

alpha 

= 

2.50000 

zB 

= 

-0.25000 

zT 

= 

0.00000 

h 

— 

0.25000 


<< Remaining elements truncated » 


Figure 8: Element parameters. 
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being quite efficient. For this manual, we will use this package for all meshing 
requirements. 

The first example is a 2 cm x 3 cm patch antenna recessed in a 5 cm x 6 
cm cavity which is centered at (<f> = 0°, z = 0 cm). The required information 
is as follows 

• Choose item 1 (Preprocessor) from main menu. 

• Enter radius of the cylinder in centimeters. 

• Enter angular and axial size of grid in degrees and centimeters. 

• Enter center of grid in degrees and centimeters. 

• Enter number of grid points in azimuthal and axial directions. 

• Enter number of identical cavities in the azimuthal and axial directions. 

• Enter number of nodes per cavity. 

• Indicate whether all surface nodes are metallic or resistive. 

• Indicate whether all surface nodes are on the substrate (0 means a 
patch is present). 

• Enter number of patches. 

• Enter row and column of lower-left hand node of the patch (see node 
section of the geometry file for this information). 

• Enter number of edges in ((f), z) directions for this patch. 

• Indicate any additional metallic nodes (-999 2 denotes no remaining 
nodes). 

• Enter number of substrate layers. 

• For each layer, enter its thickness in centimeters. 

• Enter 1 to save this geometry. 

• Enter filename for this geometry. 

The transcript for this geometry is shown on the following pages. The 
next two examples are a discrete array of these elements and a continuous 
wraparound array. The main difference is that the number of grid points is 
not the same as the number of nodes on the cavity surface. Essentially, the 
same information is entered as above with the exception that these two arrays 
have a wraparound (360°) grid. See the attached transcripts for more de- 
tails. Note for continuous wraparound cavities, if a patch crosses the branch 
cut ((f) = ±180°), the nodes along the lower edge of the patch must be 
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hand entered as shown in example 3. These nodes are obtained by running 
the preprocessor without specifying any patches and inspecting the human- 
readable geometry file. Note that for these three examples, the number of 
unknowns is dramatically different. For the single cavity of example 1, only 
541 unknowns are required. The four identical discrete cavities requires 2164 
unknowns which is of course 4 x 541. Finally, the continuous wraparound 
cavity has 12992 unknowns! 
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! i i I I ! I II I I I I I M M I I I I M I I I i I it I I I M M I i I I 

I I 

I Example 1 I 

I Single Discrete Cavity Example l 

I I 


Do you wish to run: 

1) Preprocessor 

2) FEM-BI 

3) Impedance aperture 
0) Exit 

1 

Performing initialization, please wait,... 

Enter radius of cylinder (cm) : 

15.27887 

Enter array size (phi,z) in (deg, cm): 

18.75 6 

Enter center of array (phi,z) in (deg, cm): 

0 0 

Enter number of grid points in (phi,z) direction: 
11 25 

Discretization : 

deltaPhi = 1.875 deg 

deltaZ * .25 cm 

Enter number of cavities in phi,z directions 
1 1 

Enter number of nodes per cavity (phi,z): 

11 25 


275 surface nodes have been generated. . . 

The node numbering has the following pattern 

SINGLE CAVITY: 


1 

21 

22 

23 

24 

25 

1 

16 

17 

18 

19 

20 

z 

11 

12 

13 

14 

15 

1 

6 

7 

8 

9 

10 

1 

1 

2 

• t'n Vi A 

3 

4 

5 


if for example 25 nodes were specified, 
where node 13 is the center point. 


Now specify which other nodes are either: 
metallic 
or 

resistive 

All surface nodes metal or resistive (l=yes) : 

0 

Are the nodes substrate nodes (l=yes) : 

0 

Enter number of metallic patches: 

1 

Enter row, column of lower left corner of 
the patch on uniform grid for patch: 1 
6 3 
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Enter number of unknowns (edges) for each 
dimension (phi, z) : 

4 12 

<< List of nodes constituting the patch deleted >> 

« Each of these nodes is specified by (phi,z) coord >> 

Now enter any individual nodes... 

Specify the class (metal or resistive) of all 
surface nodes which are NOT substrate nodes: 

Also indicate if 2 = metal or 3 = resistive... 

Enter node number (-999 if done) : 

-999 2 

0 resitive nodes entered. . . 

66 metal nodes entered. . . 

Enter number of substrate layers (INTEGER) : 

1 


Generating sub-surface nodes.... 
Enter layer 1 thickness (cm) : 
0.07874 

All 550 NODES have now been created 
which form 240 elements... 

Generating edges... 

1303 edges generated. . . . 


All element edges have been identified. . . 

Classifying edges. . . . 

Number of elements: 240 
Total number of nodes: 550 

Total number of edges: 1303 

Number of interior edges: 207 

Number of metal edges: 762 

Number of aperture edges: 334 

Resistive edges: 0 

Substrate edges: 334 

Number of unknowns: 541 

Determining unknown order. . . 

Do you want to save this geometry (l=yes) : 

1 

Enter EXCALIBUR filename: 
examplel . exc 

Writing binary EXCALIBUR file: 
examplel . exc 

Do you wish to run: 

1) Preprocessor 

2) FEM-BI 

3) Impedance aperture 
0) Exit 
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I i I II I I i I II I I I I i I I I I I I I I I I I ! I I M II I I II I I I 
I 

| Example 2 

| Four Discrete Cavities Example 

I 

I II I I II 1 I I I II I I I I I I I I I M I I I I 1 I i I I I I I I I I I I 

Do you wish to run: 

1) Preprocessor 

2) FEM-BI 

3) Impedance aperture 
0) Exit 


Performing initialization, please wait.... 

Enter radius of cylinder (cm) : 

15.27887 

Enter array size (phi,z) in (deg, cm): 

360 6 

Wrap-around cavity entered. . . 

Enter center of array (z) in (cm) : 

0 

Enter number of grid points in (phi,z) direction: 
192 25 

Discretization : 

deltaPhi * 1.875 deg 

deltaZ *= .25 cm 

Enter 1 if ring element, 0 if rectangular array: 

0 

Enter number of cavities in phi direction 
4 

Enter number of nodes per cavity (phi,z): 

11 25 

Enter lower left-hand (row, column) : 1 
0 43 

Enter lower left-hand (row, column) : 2 
0 91 

Enter lower left-hand (row, column) : 3 
0 139 

Enter lower left-hand (row, column) : 4 
0 187 


1100 surface nodes have been generated. . . 

The node numbering has the following pattern 

AZMUTHAL CAVITY ARRAY: 


1 19 

20 

21 

22 

23 

24 

z 13 

14 

15 

16 

17 

18 

1 7 

8 

9 

10 

11 

12 

1 1 

2 

3 

phi - 

4 

5 

6 


if for example two 3x4 cavities 
were specified 

Now specify which other nodes are either: 
metallic 


or 
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i.i mill ill hi 


resistive 


All surface nodes metal or resistive (l=yes) : 

0 

Are the nodes substrate nodes (l=yes) : 

0 

Enter number of metallic patches: 

4 

Enter row, column of lower left corner of 
the patch on uniform grid for patch: 1 
6 46 

Enter number of unknowns (edges) for each 
dimension (phi,z): 

4 12 

« Node list deleted >> 

Enter row, column of lower left corner of 
the patch on uniform grid for patch: 2 
6 94 

Enter number of unknowns (edges) for each 
dimension (phi,z): 

4 12 

« Node list deleted » 

Enter row, column of lower left corner of 
the patch on uniform grid for patch: 3 
6 142 

Enter number of unknowns (edges) for each 
dimension (phi,z): 

4 12 

<< Node list deleted >> 

Enter row, column of lower left corner of 
the patch on uniform grid for patch: 4 
6 190 

Enter number of unknowns (edges) for each 
dimension (phi,z): 

4 12 

« Node list deleted » 

Now enter any individual nodes. . . 

Specify the class (metal or resistive) of all 
surface nodes which are NOT substrate nodes: 
Also indicate if 2 - metal or 3 * resistive... 

Enter node number (-999 if done) : 

-999 2 

0 resitive nodes entered. . . 

261 metal nodes entered. . . 

Enter number of substrate layers (INTEGER) : 

1 


Generating sub-surface nodes.... 

Enter layer 1 thickness (cm) : 

0.07874 

All 2200 NODES have now been created 
which form 960 elements... 

Generating edges... 

5212 edges generated. . . . 


All element edges have been identified... 
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Classifying edges. . . . 



Number of elements: 960 
Total number of nodes: 2200 

Total number of edges: 5212 

Number of interior edges: 828 

Number of metal edges; 3048 
Number of aperture edges: 1336 

Resistive edges: 0 

Substrate edges: 1336 

Number of unknowns: 2164 

Determining unknown order... 


Do you want to save this geometry (l**yes) : 

1 

Enter EXCALIBUR filename: 
example2 . exc 


Writing binary 

EXCALIBUR file: 

example2 . exc 



Do you wish to 

run : 



1) 

Preprocessor 


2) 

FEM-BI 


3) 

Impedance aperture 


0) 

Exit 


0 


82 


I I M I I I I I I I I I I I I I I I I I I I I I I ! I I I I | I I I I I I I I 


I Example 3 

I Four Patches on Continuous 

I Wraparound Cavity Example 

I 

I I II I I II I I I I I I I M I I I I I I I I I I I II I I I I I I I I I I I 
Do you wish to run: 

1) Preprocessor 

2) FEM-BI 

3) Impedance aperture 
0) Exit 

1 

Performing initialization, please wait.... 

Enter radius of cylinder (cm) : 

15.27887 

Enter array size (phi,z) in (deg, cm): 

360 6 

Wrap-around cavity entered. . . 

Enter center of array (z) in (cm) : 

0 

Enter number of grid points in (phi,z) direction 
192 25 

Discretization : 

deltaPhi = 1.875 deg 

deltaZ - .25 cm 

Enter 1 if ring element, 0 if rectangular array: 

1 

Enter number of identical rings (>=1) : 

1 Enter number of nodes per cavity (phi,z): 

192 25 

Enter lower left-hand (row, column) : 1 
0 0 


4800 surface nodes have been generated. . . 

The node numbering has the following pattern 

SINGLE CAVITY: 

A 


1 

21 

22 

23 

24 

25 

1 

16 

17 

18 

19 

20 

z 

11 

12 

13 

14 

15 

1 

6 

7 

8 

9 

10 

1 

1 

2 

T'kVi n 

3 

4 

5 

...s 


if for example 25 nodes were specified, 
where node 13 is the center point. 


Now specify which other nodes are either: 
metallic 
or 

resistive 

All surface nodes metal or resistive (l=yes) : 

0 

Are the nodes substrate nodes (l=yes) : 

0 

Enter number of metallic patches: 
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4 

Enter row, column of lower left corner of 
the patch on uniform grid for patch: 1 
6 46 

Enter number of unknowns (edges ) for each 
dimension (phi,z): 

4 12 

<< Metallic node list deleted >> 

Enter row, column of lower left corner of 
the patch on uniform grid for patch: 2 
6 94 

Enter number of unknowns (edges) for each 
dimension (phi, z) : 

4 12 

« Metallic node list deleted >> 

Enter row, column of lower left corner of 
the patch on uniform grid for patch: 3 
6 142 

Enter number of unknowns (edges) for each 
dimension (phi, z) : 

4 12 

<< Metallic node list deleted >> 

Enter row, column of lower left corner of 
the patch on uniform grid for patch: 4 
6 190 

Enter number of unknowns (edges) for each 
dimension (phi,z): 

4 12 

This patch includes branch... 

Enter first row nodes: 

1343 

1343 176.25 -1.5 

1344 

1344 178.1249 -1.5 

1153 

1153 -180.0 -1.5 

1154 

1154 -178.125 -1.5 

1155 

<< Metallic node list deleted >> 

Now enter any individual nodes... 

Specify the class (metal or resistive) of all 
surface nodes which are NOT substrate nodes: 
Also indicate if 2 = metal or 3 s resistive... 

Enter node number (-999 if done) : 

-999 2 

0 resitive nodes entered. . . 

261 metal nodes entered. . . 

Enter number of substrate layers (INTEGER) : 

1 


Generating sub-surface nodes.,.. 

Enter layer 1 thickness (cm) : 

0.07874 

All 9600 NODES have now been created 
which form 4584 elements... 

Generating edges... 

23616 edges generated. . . . 

All element edges have been identified... 


Classifying edges.... 


Number of elements: 4608 
Total number of nodes: 9600 
Total number of edges: 23616 
Number of interior edges: 4416 

Number of metal edges: 10624 

Number of aperture edges: 8576 

Resistive edges: 0 

Substrate edges: 8576 

Number of unknowns: 12992 

Determining unknown order... 

Do you want to save this geometry (l~yes) : 

1 

Enter EXCALIBUR filename: 
example3 . exc 

Writing binary EXCALIBUR file: 
example3 . exc 

Do you wish to run: 

1) Preprocessor 

2) FEM-BI 

3) Impedance aperture 
0) Exit 
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6 Operation Modes 

EXCALIBUR has three main operation modes for FE-BI calculations (option 
2 from the main menu). They are: input impedance vs. frequency, radiation 
pattern or RCS vs. frequency and single frequency radiation and RCS pattern 
calculations. This section will describe each mode using example 1 above. 

6.1 Input Impedance - Multiple Frequencies 

The first option presented is calculation of a patch antenna’s input impedance 
at multiple frequencies. This is most useful in determining the resonant 
frequency of a patch antenna. The following information is required 

• Choose item 2 (FE-BI) from main menu. 

• Enter the stored binary geometry file. 

• Enter 1 if all elements have the same material parameters. 

• Enter complex permittivity. 

• Enter complex permeability. 

• Enter BiCG convergence tolerance, minimum and maximum number 
of iterations. 

• Enter 1 to monitor convergence. 

• Enter 1 for diagonal preconditioning and 0 for no preconditioning. 

• Enter 1 for frequency sweep of the input impedance. 

• Enter name of file to store the input impedance. 

• Enter number of probe feeds. 

• Enter location of each feed in terms of (^>=degrees, z=cm). 

• Enter which layer in which the feed is embedded. 

• Enter complex current for this feed. 

• Enter number of impedance post loads. 

• Enter frequency range (in GHz) for this sweep. 

• Return to main menu. 

The input impedance for 3.1 GHz to 3.3 GHz computed every 10 MHz is 
shown in figure 9. The transcript for computing the input impedance of an 
axially polarized 2 cm x 3 cm patch in a 5 cm x 6 cm cavity is given on the 
following page. 
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Do you wish to run: 


1) Preprocessor 

2) FEM-BI 

3) Impedance aperture 
0) Exit 

2 

Performing initialization, please wait.... 

Enter EXCALIBUR filename: 
examplel . exc 

Reading excalibur file: examplel.exe 

Material Parameter Specification. . . . 

Is the material filling constant (l=yes, 0=no) ? 

1 

Enter relative permittivity [real, imaginary ] : 

2.17 0 

Enter relative permeability [real, imaginary] : 

1 0 

Checking dimension allocations... 

Number of Boundary Integral Unknowns: 334 
Total number of UNKNOWNS: 541 

Enter tolerance, minimum and maximum iterations: 

0.01 2 500 

Do you wish to monitor convergence (0=no, 1-yes) ? 

0 

Do you want: 0 = no preconditioning, 1 = diagonal? 

1 

Do you want to compute: 

0 * RCS/Pattern, 1 * Zin, 2 - Freq. Sweep? 

1 

Enter Zin filename: 
examplel . zin 

I I I I I I I I II I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 
| PROBE FEED INFORMATION I 

I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I i I I I I I I II 

Enter number of feeds: 

1 

Enter feed location (phi,z) in (deg, cm): 

0 -0.375 

Enter layer number of feed <=: 1 

1 

Enter mag, phase of probe current (amp, deg) : 

1 0 

I I I I I I I I I I I I I I I I I I I I I II I I I I I I I I I I I I I I I I I 
| IMPEDANCE LOAD INFORMATION I 

I I I I I I I I I I I II I I I II I I I II I I I I I I I I I I I I I I I I I 

Enter number of loads: 

0 

Enter start, stop and increment freq. (GHz): 

3.1 3.3 0.01 

« Data shown in attached figure » 
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6.2 Pattern - Multiple Frequencies 

The next option presented is calculation of a radiation or RCS pattern at 
multiple frequencies. This is useful in computing the variation of gain or 
RCS with respect to frequency. Usually, a single observation angle is spec- 
ified although multiple angles are allowed. This mode permits radiation, 
bistatic and backscatter computations. In addition, the input impedance as 
a function of frequency is stored if a probe feed is used for excitation. In this 
example, we compute the backscatter at normal incidence for an E^-polarized 
plane wave as a function of frequency. The required information is 

• Choose item 2 (FE-BI) from main menu. 

• Enter the stored binary geometry file. 

• Enter 1 if all elements have the same material parameters. 

• Enter complex permittivity. 

• Enter complex permeability. 

• Enter BiCG convergence tolerance, minimum and maximum number 
of iterations. 

• Enter 1 to monitor convergence. 

• Enter 1 for diagonal preconditioning and 0 for no preconditioning. 

• Enter 2 for frequency sweep of the far-zone fields. 

• Enter name of file to store the input impedance. 

• Enter name of file to store the RCS or Gain. 

• Enter observation type (0 = backscatter). 

• Enter start, stop and increment azimuth (<f>) angles (in degrees). 

• Enter start, stop and increment elevation ( 0 ) angles (in degrees). 

• Enter polarization angle (0 = E-pol, 90 = H-pol). 

• Enter RCS filename. 

• Enter number of probe feeds. 

• Enter number of impedance post loads. 

• Enter frequency range (in GHz) for this sweep. 

• Return to main menu. 

The RCS for 3.1 GHz to 3.3 GHz computed every 10 MHz is shown in figure 
10. The transcript for backscatter calculations for backscatter calculations 
is shown on the next page. 
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Do you wish to run: 


1) Preprocessor 

2) FEM-BI 

3) Impedance aperture 

0) Exit 

2 

Performing initialization, please wait.... 

Enter EXCALIBUR filename: 
examplel . exc 

Reading excalibur file: examplel.exe 

Material Parameter Specification. . . . 

Is the material filling constant (l«yes, 0=no) ? 

1 

Enter relative permittivity [real , imaginary ] : 

2.17 0 

Enter relative permeability [real, imaginary ] : 

1 0 

Checking dimension allocations... 

Number of Boundary Integral Unknowns: 334 
Total number of UNKNOWNS: 541 
Enter tolerance, minimum and maximum iterations: 

0.01 2 500 

Do you wish to monitor convergence (0=no, 1-yes) ? 

0 

Do you want: 0 = no preconditioning, 1 *= diagonal? 

1 

Do you want to compute: 

0 - RCS/Pattern, 1 * Zin, 2 = Freq. Sweep? 

2 

Enter Zin filename: 
test . zin 

Enter Frequency sweep data filename: 
examplel . f sw 
Enter observation type: 

0) Backscatter 

1) Bistatic 

2) Radiation 


Enter start, stop, and increment azmuth angles [deg]: 

0 0 1 

Enter start , stop, and increment elevation angles [deg] 
90 90 1 

Enter polarization angle [0 <■ alpha <* 90 deg] : 
Ephi=0: alpha = 0 

Etheta=0: alpha *= 90 

0 

Enter ouput RCS filename [<= 40 characters] : 
junk 

I I I i I I I I I I I I I I 1 I I I I I I I I I I II I I I I I I I I I I I I I 

| PROBE FEED INFORMATION | 

I I I I I I I I II I I I I I I I II I I I I I I I I I I I I I I I I I J I I [ 

Enter number of feeds: 

0 

I I I I I II I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

| IMPEDANCE LOAD INFORMATION I 

I I II I I I I I II I I I I I I I I I I I I I I I I I II I I I I I I I I I I 

Enter number of loads: 

0 

Enter start, stop and increment freq. (GHz): 

3.1 3.3 0.01 


6.3 Pattern - Single Frequency 

The final operation mode is radiation and RCS pattern calculations at a 
single frequency. Of course, for the case of an antenna, the input impedance 
is also computed. This mode is generally used for multiple incident and 
observation angle applications. The only difference between this mode and 
the previous two modes is the option to save the admittance matrix in binary 
format. Since this matrix is excitation independent (though still frequency 
dependent), it can save time for multiple runs of a large geometry if this 
matrix is stored and then read in for each run. The required inputs are 

• Choose item 2 (FE-BI) from main menu. 

• Enter the stored binary geometry file. 

• Enter 1 if all elements have the same material parameters. 

• Enter complex permittivity. 

• Enter complex permeability. 

• Enter BiCG convergence tolerance, minimum and maximum number 
of iterations. 

• Enter 

• Enter 

• Enter 

• Enter 
ing). 

• Enter 

• Enter 

• Enter 

• Enter 

• Enter 

• Enter 

• Enter 

• Enter 

• Enter 

• Enter 

• Return to main menu. 

This option is illustrated in the next session transcript for bistatic scattering 
and radiation pattern calculations. 


1 to monitor convergence. 

1 for diagonal preconditioning and 0 for no preconditioning. 

0 for single frequency operation. 

matrix storage/read option( 1 to store, 2 to read, 0 to do noth- 

six (6) character filename for matrix storage/read, 
observation type (1 = bistatic), 
incidence angle (<f>,0). 

start, stop and increment azimuth ( <f > ) angles (in degrees), 
start, stop and increment elevation (0) angles (in degrees), 
polarization angle (0 = E-pol, 90 = H-pol). 

RCS filename, 
number of probe feeds, 
number of impedance post loads, 
frequency (in GHz). 
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I PIP 



Do you wish to run: 


1) Preprocessor 

2) FEM-BI 

3) Impedance aperture 

0) Exit 

2 

Performing initialization, please wait.... 

Enter EXCALIBUR filename: 
examplel . exc 

Reading excalibur file: examplel.exe 
Material Parameter Specification. . . . 

Is the material filling constant (l=yes , 0=no) ? 

1 

Enter relative permittivity [ real, imaginary ] : 

2.17 0 

Enter relative permeability [ real, imaginary ] : 

1 0 

Checking dimension allocations... 

Number of Boundary Integral Unknowns: 334 
Total number of UNKNOWNS: 541 
Enter tolerance, minimum and maximum iterations: 

0.01 2 500 

Do you wish to monitor convergence (O^no, 1-yes) ? 

0 

Do you want: 0 = no preconditioning, 1 = diagonal? 

1 

Do you want to compute: 

0 * RCS/Pattern, 1 - Zin, 2 = Freq. Sweep? 

0 

Save the BI arrays in a file (0=no, l=yes, 2=read) ? 

1 

All BI files will have the following form: 

<f ilename> . guu, <f ilename> . guv, etc . 

Enter BI binary filename (must be six (6) char): 
exampl 

Enter observation type: 

0) Backscatter 

1) Bistatic 

2) Radiation 

1 

Enter incident phi, theta [deg]: 

0 90 

Enter start, stop, and increment azmuth angles [deg]: 
-180 180 1 

Enter start , stop, and increment elevation angles [deg] 
90 90 1 

Enter polarization angle [0 <= alpha <= 90 deg] : 
Ephi=0: alpha = 0 

Etheta=0: alpha = 90 

0 

Enter ouput RCS filename [<*= 4 0 characters] : 
exampl . bi 

Enter frequency (GHz) : 

3 . 3 

I I I I I I I I I I I II I I I I I I I I I I I I I II I II I I I I I I I I I 

| PROBE FEED INFORMATION I 

I I I I I I I II I I II I I I I I I I I I I I I I I I I I I I I I I I I I I I 
Enter number of feeds: 

1 

Enter feed location (phi,z) in (deg, cm): 

0 -0.375 

Enter layer number of feed <=: 1 

1 

Enter mag, phase of probe current (amp, deg): 

1 0 
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I I III I ! I M I i II II I I II I I I II I I I I II ! I I I I I I I I 

i IMPEDANCE LOAD INFORMATION I 

I M M I I II ! I I I I I (I I I I I I I I 11 I ! I II [ i I ! II I I ! 

Enter number of loads: 

0 

Writing BI binary files... 

Generate bistatic excitation vector... 

Solve it . . . 

« Pattern deleted. » 

Run anouther excitation? (l=yes, 0=no) 

1 

Enter observation type: 

0) Backscatter 

1) Bistatic 

2) Radiation 

2 

Enter start, stop, and increment azmuth angles [deg]: 
-180 180 1 

Enter start , stop, and increment elevation angles [deg]: 
90 90 1 

Enter ouput RCS filename [<= 40 characters] : 
exampl . pat 

Enter normalized pattern filename: 
exampl . norm 
Solve it,., 

« Pattern deleted. » 

Run anouther excitation? (l~yes,0=no) 

0 

Do you wish to run: 

1) Preprocessor 

2) FEM-BI 

3) Impedance aperture 
0) Exit 
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The bistatic patterns taken at 3.3 GHz for the three geometry configu- 
rations presented in this manual are shown in figure 11 The corresponding 
antenna pattern comparisonat 3.3 GHz is shown in figure 12 


7 Concluding Remarks 

This user manual presented some basic operation information for the FE- 
BI code, EXCALIBUR. This presentation was only meant to get an initial 
user started. As one becomes experienced with the code, additional features 
such as 2-D patch array modeling, multiple feed arrays and use of lumped 
impedance loads may prove useful. Indeed, an experienced user will find 
that custom features may readily be added to EXCALIBUR. For example, 
currently EXCALIBUR allows entry of material parameters either for the 
entire substrate, each layer of the substrate or on an element-by-element 
basis. This subroutine material in file excalibur.f may readily be modified by 
the user to input a custom inhomogeneous substrate. 

The code is fairly “dummy proof”. If the user enters data which is not 
expected by EXCALIBUR such as a character when an integer is expected or 
an angle greater than 360°, the code will prompt the user to re-enter the re- 
quested data. Additionally, the storage allocation parameters in excalibur.inc 
must be set by the user prior to compilation. If a particular parameter is 
too small for a given run, the code will halt and suggest a new value for the 
offensive parameter. The user must reset that parameter, recompile and run 
the code again. The program also estimates the amount of RAM required 
at the start of a run. This estimate is based on the storage required by the 
arrays in excalibur.inc plus some scratch arrays. Each complex number is 
assumed to require eight bytes and each integer and real number require four 
bytes. The user should consider this estimate to be a slightly lower than the 
actual consumed memory. 
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Frequency (f) [GHz] 


Figure 10: RCS vs. frequency for a patch antenna which is 2 cm x 3 cm in 
a 5 cm x 6 cm cavity. 






Angle (<J>) [deg] 


Figure 11: E-polarized bistatic patterns for a 2 cm x 3 cm embedded in the 
three different geometries presented in this manual. The incidence angle is 
(0°,90°) and observation is in the 6 — 90° plane. 
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Angle (<|>) [deg] 


Figure 12: Axially polarized H-plane patterns for a 2 cm x 3 cm embedded 
in the three different geometries presented in this manual. 




Appendix 1 


A Note ox Applying New Vector ABCs to Patch 
Antennas on a Circular Cylinder 


Leo C. Kempel and John L. Volakis 
Radiation Laboratory 

Department of Electrical Engineering and Computer Science 
University of Michigan, Ann Arbor, MI 48109-2122 

Recently, we developed a Finite Element-Boundary Integral (FE-BI) method which 
may be used to accurately model the scattering and radiation by conformal patch 
antennas. The use of a boundary integral to provide an exact mesh closure condition 
results in a rigorous formulation. When coupled with a Biconjugate Gradient-Fast 
Fourier Transform (BiCG-FFT) solver, the resulting computer code possesses both 
low memory and low computational demand. However, boundary integral closure 
involves a dyadic Green’s function which necessarily limits its extension to coated 
geometries. Although it is possible to construct a suitable Green’s function for mul- 
tilayered superstates, the resulting formulation is both complex and bulky. Fur- 
thermore, this approach cannot be extended in a rigorous manner to more general 
doubly curved structures. 

Another hybrid FE method which has been used for very large 3-D scattering 
calculations at the University of Michigan is the finite element-absorbing boundary 
condition (FE-ABC) method. ABCs provide an approximate mesh closure which 
maintains a highly sparse matrix for the whole system. We propose utilizing a 
new conformal ABC (A. Chatterjee and J.L. Volakis, Microwave and Optical Tech. 
Letters , 6, No. 16, pp. 886-889, Dec. 20 1993) which allows the closure surface to 
be brought quite close to the cavity aperture. As a result, the number of unknowns 
required is significantly less than the number required if a traditional spherical 
closure surface/boundary condition is used. 

We have applied the second order ABC to radiation and scattering by cavity- 
backed patch antennas which are embedded in a metallic cylinder. Although ex- 
tensive testing has not as of yet been accomplished, the preliminary results are 
encouraging. We looked at a 2 cm x 3 cm patch antenna which was printed on a 

5 cm x 6 cm x 0.07874 cm substrate which had a dielectric constant of e r = 2.17. 
The second order ABC was place 0.5A from the cavity aperture as shown in figure 
1. Two comparisons with the FE-BI formulation were performed. The first involved 
bistatic scattering with normal incidence (<£, = O°,0; = 90°) and observation in the 

6 = 90° plane for an ^-polarized plane wave which is shown in figure 2. A radiation 
pattern for the same antenna is shown in figure 3 where the resonant frequency is 
3.3 GHz. Since this is an axially polarized element, figure 3 is an H-plane pattern. 

The scattering and radiation example given herein involves a conformal antenna 
element without a covering superstrate. The FE-ABC formulation may be readily 
modified to account for such a covering which may consist of multiple inhomogeneous 
layers as shown in figure 1. A future report will develop the formulation for both 
covered and uncovered elements and it will provide suitable validation of the method. 
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In addition, we shall investigate whether the near-field performance of the FF-ABC 
method is sufficient for input and mutual impedance calculations. Having shown 
that the FE-ABC approach can indeed be used for conformal antennas on singly 
curved surfaces, we may in confidence develop a similar formulation for antennas 
mounted on doubly curved surfaces. 



Figure 1: Coated cavity-backed patch antenna with ABC mesh termination. 


100 


I liiiliiinii 




Figure 2: Comparison of 2 nd order conformal ABC mesh closure condition 
with the exact boundary integral condition for bistatic scattering. The patch 
was 2 cm x 3 cm printed on a 5 cm x 6 cm x 0.07874 cm dielectric substrate 
(e r = 2.17). The ABC was placed 0.5A from the aperture and the operating 
frequency was 3.0 GHz. The ^-polarized plane wave was incident from (</>, = 
O°,0, = 90°) and observation is in the 0 = 90° plane. 



Angle (<}>) [deg] 

Figure 3: Comparison of 2 nd order conformal ABC mesh closure condition 
with the exact boundary integral condition for antenna pattern calculations. 
The geometry is identical to the one used in figure 2. The operating frequency 
was 3.3 GHz. The feed point was (<f> s = 0 °,z s = —0.375 cm) which results in 
axial polarization and observation is made in the H-plane. 
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APPENDIX 2 


FE-BI implementation for coated planar platforms 

Consider the cavity recessed in a coated ground plane as shown in Fig 1. The 
coating is assumed to be of thickness d and of relative permittivity e T . We 
will denote the free space region above the coating {z > d ) as region 1, the 
region inside the coating (0 < z < d) as region 2 and that inside the cavity 
(— c < 2 < 0) as region 3. We will assume that the cavity is filled with 
an inhomogeneous material having a relative permittivity f,. c (r) and relative 
permeability /z rc (r). In accordance with the equivalence principle, the fields 
in regions 2 and 3 can be decoupled by closing the aperture with a perfect 
conductor and introduce the equivalent magnetic current 


M = E x z 


( 1 ) 


over the extent of the aperture, where E is the electric field at the aperture ( z = 
0). The magnetic field is then due to the radiation caused by the equivalent 
current M residing on the ground plane and is given by 


H(r) 


'(r) + H r (r) — jkY fl + -^W 


( 2 ) 


where F is the magnetic vector potential and is given in terms of the dyadic 
Green’s function of the coating as 


F = y^G(r,r')-M(r')dS' 


(3) 


and H‘ is the incident field, H r is that reflected by the coated ground plane 
without the aperture, S denotes the planar surface area of the aperture. If 
the conducting surfaces are restricted to horizontal planes, as is the situation 
here, we can find an equivalent scalar Green’s function Gv such that [1] 

V • G = V GV(r,r') (4) 


where 


Gv 


u dfff 
kjk x dz 



(5) 
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More specifically, the function P defined in [2] is identically zero. Now consider 

F s as 


Fs 


—jkY 


r + -w 


F 


(6) 


Making use of (4) in (6) we can write Fs as 

F s = -jkY JJ G • M(r ')dS' + V JJ^ V' • M(r')Gv(r, r ’)dS' (7) 


The required components of the Green’s function in (7) are given by 


AO O fOO 

G xx = / G xx t- ik ^ s - d h ik ^ x - x 'h ik ^ y - y 'Uk x dk y 

J — oo J — oo 

( 8 ) 

r oo r oo 

G yx = / / G yx e-^ {z - d) e ikl{x - x ">e ik ^ y - y,) dk x dk y 

J — OO */ — OO 

( 9 ) 

/“OO /“OO 

G w = / / G yv e- A(2 - J) e A(l - 1 >eJ'M»-v )<ffc x d£ y 

*/ — oo */ — oo 

(10) 

A OO /“OO 

G xy = / / Gl y e~ jk ^ z - d) e jk ^ x - x,) e jk ^ y - y,) dk x dk y 

J — OO J —oo 

(11) 

where 

<?« = / (8 

(12) 


g„ = / {zxXhf) 

(13) 


<e 

ii 

P 

H 



= -fc x fc v (fc 2 T 2 + Cr^Tx) / (87T 2 fcofcx^ 2 ) 

(14) 

with r x and T 2 

given by 



e r fc 2 sin(fcxd) — jk x cos(fcxd) 
1 e r fc 2 cos(fcxd) + j ki sin(fcxd) 

(15) 


fc 2 cos ( &x d) + j k-i sin( &x d) 
fc 2 sin(&id) — jk\ cos(kid) 

(16) 

Also, 

Gf s = r f°° G~f x e jkl(x - x '^e jk ^ y - y,) dk x dk y 

J -oo J — OO 

(17) 
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/■ oo r oc 

G'J X = / G^e^-^cMy-v'UkJky 

J — oo J — oc 

(18) 

with 

a\ = jlj [C «*(!,*) + Dsin(V-)] 

(19) 

where 

ky [&i cos(k\d) 4- jt r k 2 siu(A.id)] 
k i [e r k 2 cos (kid) -|- jki sin( A:! <f)] 

(20) 


D = — 
k i 

(21) 

and 

G*? r = [E sin(fciz) + Fcos(fciz)] 

(22) 

where 

k x [ki sin(fcid) — jk 2 cos(&id)] 
k 0 Z 0 [fci cos(&id) + jk 2 sin(fcjd)] 

(23) 


II 

O 

(24) 

with k 2 = 

e T k,l — P 2 , k% = k% — p 2 and P 2 = k 2 + k 2 . After 

weighting (7) 


becomes 


F s = —jkY Jj M(r) ■ JJ s G M(r ')dS dS' 

Us V ' Us V ' M ^ Gv ^ T ') dS dS ' ( 25 ) 

The fields in region 3 (the cavity region) are formulated employing the finite 
element method [3]. 

Evaluation of the Sommerfeld integrals using the Chebyshev de- 
composition method 

This is a new, fast and efficient method [4] of evaluating Sommerfeld integrals 
encountered, when we use a spectral-domain transformation to calculate fields 
of a dipole source in a homogeneous space or in a layered medium with planar 
boundaries. The integration is based on the decomposition of the integrand 
into Chebyshev polynomials, the oscillatory nature of which facilitates rapid 
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integration of the fast oscillating integrands. In (S)- (18) we encounter Fourier 
transform pairs given by 



1 f 

/(*«, k y ) = — / / f(x, y)e~* kxZ ~ }kyV dxdy 

Z'K J J — oo 

(20) 


1 f°° t 

f{x,y) = — / / f(k x , k y )e^ x+ ^ dk T dk y 

Ztt J J — oo 

(27) 

Introducing the Bessel function as 



J 0 (k pP ) = r e^ kppco ^ a ~^d<f> 

zn Jo 

(28) 

we have 

foo 

f(k p ) = / J 0 (k p p)f(p)pdp 

Jo 

(29) 


foo 

fp= Jo{k p p)f(kp)k p dk p 

Jo 

(30) 


(29) and (30) is the Fourier- Bessel or Hankel integral transform pair. Inverse 
Hankel tranforms such as (30) are the generic form of Sommerfeld integrals 
which we consider for the demonstration of this method. Transverse derivatives 
of f(p) can be expressed as : 

d f f°° 

— = - cos <t> h J l {k p p)f{k p )k 2 p dk p (31) 

Hence the most general Sommerfeld integral is written as 

SI = f°° J m (k pP )h(k p ) dk p (32) 

J o 

where J m is the Bessel function of the first kind of order m. Assuming that 
h{k p ) has a pole at k p = kp (0 < kp < y/e^k 0 ) 

fy/trh o 

SI = / (J m (k pP )h(k p )-F(k p ))dk p 

J o 

1 F(k p ) dk p 

0 

[OO 

+ / Jm{k p p)h(k p ) dkp (33) 
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with 


F(k P ) = 


A 

k p k p 


(3-1) 


where A is the residue of J m (kpp)h(k p ) at k p = kp. The second integral in (33) 
is solved analytically. Next we make the behaviour of the integrands in (33) 
smoother by applying a change of variables, for the first and third integrands 
respectively 

k p = y/T r k Q sin i (35) 

and 

kp = y/^rkf) cosh £ (36) 

Transforming the domain of integration to the domain of definition of the 
Chebyshev polynomials the first integral of (33) becomes 


u' 4hU 1+ ^/' 

4 J - 1 y/l ~z\ 4 J ~ l \,A- 


A 


(37) 


with 


= Re , Im (( J m (k pP )h(kp ) - F(k,))y/l - kjy/l - Z 2 ) (38) 

kp = y/T r k 0 sin + 1)^ (39) 

Similarly transforming the domain of integration for the third integral in (33) 
and truncating the upper limit at a and b for the real and imaginary parts 
respectively we have 


A [ l hMfsMj , ;t> f l UWJi 

2L-jTzT dZ2+ ^J-,-Jrz 


.b r l /tfeO/sfe) 




v/i 


dz?t 


(40) 


with 


fsA 2 ) = J m (y^ko cosh (^y(z + 1)^ pj 


fs.ei 2 ) = R e 


, Im ^h, ^ye^ko cosh ^ 


ij j y/l - z 2 

(41) 


(42) 
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(43) 


Now f\(z) to fe(z) an 1 fi i ted with Ohebyshrv polynomials so that 

A r 


/(*) 




fc=i 


1 

-Ci 

2 


with 


1 k~l 

The orthogonal property of these polynomials is defined as 


L 


1 T^T^z) 




o (i ^ j) 

* = { f (i=i*0) 

7T (i = j = 0) 


Applying the orthogonal property we evaluate the SI as 

m 7T7T x 7T ( a v b \ 

SI ~ 2 4 ^ C1,1 +;C2 ’^ + 2 \ 2 ^ C3,C5 ^ + - 7 2^ C4,C6 7 

/ c T k 0 kp 


+A\n 


V k P 


— jit A 


with 


/(ci ' c>)= £ ■ 2 h 




(44) 


(45) 


(46) 

(47) 


Using the Chebyshev decomposition method, a typical Sommerfeld integral 
encountered while analysing a horizontal electric dipole in the x direction 
above a ground with relative pemittivity e r and conductivity <r, is evaluated. 
The integral is expressed as 

/■oo / g — uo|*+fc| \ 

A ’ =2 l (-^nr j kfikf (48) 

where 

P = \A 2 + y 2 u 0 - - k 2 (49) 

u = M - iW N* = e ,-j— (50) 

o;eo 

and A: is the free space propagation constant. For this evaluation z (the ob- 
server height) to be equal to the dipole height h (A/50). The results of this 
evaluation is shown in Fig 2. 
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Region 1 



Region 2 Region 3 


Figure 1: Geometry of a groove in a coated plane 
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(V)sqv 


Variation of the Sommerfeld potential A . of a HED 

X 



Figure 2: A typical Sommerfeld integral evaluation 
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